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Synopsis 


The dynamic characteristics of the rotor blades play a significant role in the overall 
performance and stability of a helicopter. There has been a continued effort to develop 
a mechanically simple yet efficient rotor blade and hub configuration. Early rotor 
blades were provided with flap (out-of-plane) and lag (in-plane) hinges at the root of 
the blade. These rotor systems having complicated construction are usually referred 
to as articulated rotors. With the advancement in fibre-reinforced composite material 
technology', increasing emphasis has been placed on the development of hingeless and 
bearingless rotor systems. Both articulated and hingeless rotors are provided with 
external hydraulic dampers to increase damping in the lag mode and thereby to avoid 
aeroelastic and/or aeromechanical instabilities. In a bearingless rotor, external damper 
is eliminated by incorporating a specialised elastomer with high loss factor. 

For the last two decades, considerable effort has been put in by industry and R &: D 
establishments to develop a bearingless main rotor. Though a bearingless rotor is me- 
chanically simple, its dynamic analysis becomes very complicated because of multiple 
load paths, presence of non-linear elastomeric damper and existence of kinematic con- 
straint at the pitch link. With the consolidation of rotor blade structural dynamic 
and aeroelastic formulation applicable for articulated and hingeless rotor blades, aca- 
demic research has now focussed on the fundamental understanding of the effects of 
non-linearities of elastomer on the dynamics of bearingless rotor blades. 



In general, the elastomer used in a bearingless rotor exhibits highly non-linear char- 
acteristics both in stiflfness and damping with respect to the amplitude of deformation. 
It is expected that due to the presence of the elastomer, the lag and flap frequencies 
should dependent on the amplitude of motion which is a typical feature of a non-linear 
system. But so far no attempt has been made to study the effect of amplitude on the 
frequencies of a bearingless blade. Although several researchers have proposed different 
non-linear models for the elastomer, while performing the stability analysis, only lin- 
earised equations of the blade and elastomer are solved. In addition, all the aeroelastic 
response studies of the rotor blade have focussed on the steady state response. No 
information is available on the trasnsient response characteristics of a bearingless rotor 
with the inclusion of kinematic constraint and multiple load paths. 

The major objectives of the present study are as follows. 

(i) Formulate a simple non-linear model to capture the non-linear stiffness and damping 
characteristics of an elastomer. 

(ii) Study the influence of non-linear stiffness of the elastomer on the free vibration 
characteristics of uncoupled flap and lag dynamics of the blade. 

(iii) Examine the amplitude dependent stability of an idealised bearingless rotor system 
under ground resonance. 

and 

(iv) Analyse the transient response of an isolated bearingless rotor blade undergoing 
coupled flap-lag-torsional deformations in hovering condition. The transient response 
results are compared with those of a hingeless blade with identical blade properties. 

A brief outline of the chapterwise contents of the thesis is summarised below^ 

The experimental data on the non-linear stiffness and damping properties of an 
elastomer show very weak dependence on frequency in the range of interest for rotor 
blade analysis. In the light of this observation, the non-linear characteristics of the 
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elastomer is modelled in Chapter 2 by a parallal combination of a non-linear spring, 
a Coulomb damper and a hysteretic damper. The model is kept in a simple form so 
that it cam be easily integrated with the blade model. The pairameters of the analytical 
model are obtained by correlating with the experimental data available in the open 
literature. 

In Chapter 3, both linear and non-linear free vibrations of an idealised bearingless 
rotating blade undergoing uncoupled flap and lag bending have been studied. First 
a linear problem hzis been solved using two different solution techniques, one based 
on power series expansion and the other based on the Rayleigh-Ritz method. The 
natural frequencies are obtained for different values of spring stiffnesses and elastomer 
locations. In the non-linear analysis, a numerical-perturbation technique based on 
multiple-time-scale is formulated to determine the frequency-amplitude relationship 
for the rotating blade. 

In Chapter 4, an amplitude dependent stability analysis has been carried out for a 
coupled rotor/fuselage system under ground resonance condition. The stability analysis 
is carried out for different locations and amplitudes of motion of the elastomer. 

In Chapter 5, a suitable numerical technique has been formulated to study the 
transient response of an isolated bearingless rotor blade to a step control pitch input. 
The blade is assumed to undergo coupled flap-lag-torsional deformation in hovering 
condition. A time varying inflow model (based on the dynamic inflow model) has been 
used for the analysis. The expressions for aerodynamic loads are treated by an implicit 
formulation. The study is carried out for different locations of the elastomer and pitch 
link. The transient response characteristics of the bearingless rotor blade is compared 
with that of the hingeless blade. The comparison is made to highlight the difference in 
their response behaviours. A quantitative measure of the aeroelastic couplings (pitch- 
flap, pitch-lag), based on the steady state response has been proposed and an attempt 
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is made to relate the couplings and the response characteristics of the blade. 

An analysis has been carried out, in Chapter 6, to investigate the phenomenon of 
limit cycle oscillation and elastomer model. Finally, in Chapter 7, the overall con- 
clusions are included. The major conclusions of the present thesis are summarised 
below. 

(i) Using an idealised model for a bearingless rotor, both linear and non-linear free 
vibrations of the rotating blade have been studied. Two different solution techniques 
are used for the linear analysis. The natural frequencies obtained by both techniques 
are in excellent agreement. The results of the non-linear free vibration analysis of 
the rotor blade indicate that upto a fairly high-value of amplitude, a seventh order 
approximation of the non-linear spring is sufficient to correctly predict the frequency- 
amplitude relationship. 

(ii) It is revealed from the ground resonance stability analysis that the elastomer am- 
plitude does not significantly alter the stability of the system in the region of maximum 
instabilities. However, it exhibits appreciable influence in other regions. The elastomer 
locations seem to have considerable influence on the stability of the system. The results 
indicate that for a given location of the elastomer, there is an optimum location for 
the torque tube attachment. 

(iii) The transient response analysis to a step control pitch input is performed for both 
hingeless and bearingless rotor blade configurations having identical blade properties. 

It is observed that there is a qualitative difference in the nature of the response as well 
as in the dynamic overshoot of the rotor thrust for the two types of blades. 

(iv) It is noticed that the transient response characteristics of a bearingless blade de- 
pends on the pitch link locations. In particular, it is observed that the blade response 
with trailing edge pitch link configuration is significantly different from that of a lead- 
ing edge configuration. For certain trailing edge configurations, the blade exhibits a 
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divergence transient response indicating instability. The reason for these significant dif- 
ferences in the response behaviour can be attributed to the changes in the aeroelastic 
couplings and the constraint loads due to the pitch link. The aeroelastic coupling mea- 
sures are observed to be dependent on the magnitude of the control input, exhibiting 
a non-linear relationship. 

(v) The results of the transient response characteristics with two different elastomer 
models (for the same experimental data) indicate that the phenomenon of limit cycle 
oscillation is highly dependent on the elastomer model. Therefore, sufficient care must 
be excercised in modelling the elastomer. 
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Chapter 1 


INTRODUCTION 

1.1 Introduction 

The main rotor system of a helicopter is the most important component of the vehicle. 
The rotor system has to fulfil multiple roles of generating lift, propulsive and control 
forces. Therefore, the dynamic characteristics of rotor blades play a significant role 
in the overall performance and stability of the vehicle. There has been a continued 
effort to develop a rotor blade and hub configuration which is mechanically simple yet 
efficient. During operation, the blades experience large bending and centrifugal loads. 
In order to relieve the root bending moments experienced by the blades, early rotor 
blades were provided with flap (out of plane bending) and lag (in-plane bending) hinges 
at the root of the blade. In addition, a pitch control bearing was provided to control 
the pitch angle of the blade. Such rotor systems are usually referred to as articulated 
rotors. A schematic diagram of an articulated rotor system is shown in Fig. 1.1a. While 
the provision of hinges represents an efficient engineering solution to the problem of 
alleviating high bending loads at the root, there are several disadvantages of this type 
of arrangement. The large number of moving parts leads to a mechanically complex 
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rotor hub system, accompanied by the associated wear out problem requiring frequent 
maintenance and replacement of parts. The presence of hinges prevents generation of 
large control moments and thereby restricts the flight envelope and the range for C.G. 
travel in the helicopter. The mechanical complexity of an articulated rotor hub can be 
seen in Fig. 1.1b (taken from Ref. [1]). 

With the advancement in fibre-reinforced composite material technology, increasing 
emphasis has been placed on the development of hingeless rotor systems. Figures 1.2a 
and 1.2b show respectively a schematic and an actual rotor hub configuration of a 
hingeless rotor. The construction of these rotors is relatively simple because of the 
absence of flap and lag hinges; but a pitch bearing is still provided for blade pitch 
control. Due to the absence of hinges, large control moments can be generated which 
in turn provide favourable control characteristics of the vehicle. Both articulated and 
hingeless rotors are provided with external hydraulic dampers to increase damping in 
the lag mode and thereby to avoid aeroelastic and/or aeromechanical instabilities. 

The development of a bearingless rotor system aims to eliminate both the pitch 
bearing as well as the external damper by incorporating a specialised elastomer with 
high loss factor. In this rotor system (Fig. 1.3a), the blade is attached to the hub 
through a flexible structural element called flexbeam. The flexbeam is designed to 
provide the required stiffnesses in the flap and lag bending deformations of the blade, 
but it is highly flexible in torsion. Surrounding the flexbeam, there is a stiff cuff denoted 
as torque tube which is attached to the blade-flexbeam junction at the outboard end 
and to a pitch link at the inboard end as shown in Fig. 1.3a. The pitch control of 
the blade is achieved by rotating the torque tube through up/down movement of the 
point P which in turn twists the flexbeam. An elastomer is placed between the torque 
tube and the flexbeam to provide adequate lag damping. It also serves the purpose 
of a spacer between the torque tube and the flexbeam. Though a bearingless rotor is 
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mechanically simple (Fig. 1.3b), its dynamic analysis [2] becomes very complicated 
due to the presence of multiple load paths, non-linear elastomeric damper and existence 
of a kinematic constraint at the pitch link. 

Over the years, considerable effort has been put in by the helicopter industries 
in developing soft-in-plane bearingless main rotors with acceptable performance and 
stability characteristics. Huber [3] presented a comprehensive review of the bearingless 
main rotor programme pursued at different helicopter industries in the world. A brief 
history on the development of bearingless rotors is provided in Refs. [4, 5]. Lockheed 
developed the first protot}T)e of a bearingless rotor in mid 1960’s having “matched 
stiffness” (equal stiffness in both in-plane and out of plane bending). During test 
flight, the prototype encountered an air-resonance instability which was eliminated 
by increasing the fuselage inertia. During late 1970’s and early 1980’s, intense activity 
continued in developing the prototype of a Bearingless Main Rotor (BMR). The Boeing 
Helicopter Company developed a BMR with a high fundamental flap frequency and 
tested it on an MBB BO-105 Helicopter. It was observed that the landing gear of 
BO-105 had to be stiffened in order to improve the ground resonance stability. Another 
design was proposed by Bell Helicopter with a four bladed M680 rotor. An elastomeric 
damper was provided to augment edgewise damping. In 1990, the Boeing-Sikorsky 
Company developed a five-bladed BMR (for the new helicopter known as RAH-66 
Comanchi) with elastomeric dampers. This vehicle is currently undergoing tests. In 
Europe, Eurocopter recently built and flight tested EC135 helicopter with a BMR 
[6]. For the analysis and design of bearingless rotor blades, most of the rotorcraft 
manufacturers have developed their own in-house methodologies and computer codes. 
Hence, very little information is available about them in the open literature. 

The structural dynamic and aeroelastic/aeromechanical analyses of bearingless ro- 
tors require the formulation of suitable analytical models and solution techniques. 
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These mathematical models must include the following important features: 

i) non-linear stiffness and damping characteristics of the elastomer, 

ii) kinematic constraint at the pitch link 
and iii) multiple load paths. 

Several research studies have been reported in the open literature on modelling and 
analysis of bearingless rotor blades. A survey of the available literature is provided in 
the following section. 

1.2 Literature Review 

The literature review is divided into four subsections with each one addressing a par- 
ticular topic relevant to the study of bearingless rotor blades. 

1.2.1 Elastomeric Damper 

Elastomeric materials are becoming increasingly popular for vibration attenuation ap- 
plications due to the tremendous advantages they offer over mechanical dampers. In 
the case of bearingless rotors, conventional hydraulic dampers are being replaced by 
specialised elastomeric material with high loss factor. It is well known that under 
dynamic conditions, elastomeric materials exhibit visco-elastic behaviour dissipating 
energy through hysteresis. Besides exhibiting frequency and temperature dependence, 
the characteristics of an elastomeric material are also non-linear with respect to the am- 
plitude of motion. Consequently, it has been a challenging task to develop an analytical 
model that can accurately represent the behaviour of the elastomeric dampers. 

For visco-elastic materials, the stress-strain relations depend on the time history. 
Classical linear constitutive relations for visco-elastic materials are formulated either 
in the time domain or in the frequency domain [7]. In the frequency domain approach, 
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the coefficients of complex moduli are obtained by fitting a curve to the steady-state 
experimental data obtained over a limited frequency range of interest. 

In the time domain approach, the constitutive relationship can be formulated in 
two different ways [7]; (i) differential equation form involving derivative operators, 
and (ii) convolution integral approach. In the differential equation form, higher order 
time derivatives of stress and strain are required to represent the dynamic properties of 
material over a broad range of frequency. On the other hand, the integral formulation 
avoids the presence of higher order derivatives and is very useful for creep/relaxation 
studies of visco-elastic materials. However, in structural dynamic problems represented 
by a set of differential equations, the incorporation of integral representation of visco- 
elastic material behaviour leads to integro-differential equations which are cumbersome 
to solve. 

The method of irreversible thermodynamics has been used by Lesieutre and Mingori 
[8] to develop coupled material constitutive relations. These equations are used to 
model frequency dependent material damping in finite element domain. Later, a time 
domain model for linear visco-elastic materials has been proposed by Lesieutre and 
Bianchini [9] using Anelastic Displacement Field (ADF). The Golla-Hughes-Mctavish 
mini-oscillator model [10] is a similar type of visco-elastic formulation where the higher 
order derivatives are replaced by internal dissipation coordinates. 

Instead of integer derivative operators in the constitutive equation, introduction of 
fractional derivative operators can reduce the number of terms and parameters required 
to model visco-elastic materials. This fractional derivative approach [11, 12], being 
a time domain approach, is suitable for transient analysis. Bagley and Calico [13] 
and Fenander [14] developed a finite element formulation using fractional derivative 
approach. In a later study [15], this formulation was applied to analyse the transient 
response of a damped beam. Recently, Enelund and Josefson [16] used the fractional 
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derivative approach to capture the weak frequency dependence of the loss factor of an 
elastomer. 

Apart from the above described linear models, many theories [17]- [19] have been 
proposed to model the non-linear behaviour of visco-elastic materials using hereditary 
integral representation. Glockner and Szyszkowski [20, 21] proposed a semi-empirical 
constitutive model to predict creep, strain softening and relaxation behaviour. These 
non-linear models, being in the integral form, cannot be conveniently incorporated into 
the structural dynamics analysis. 

Due to the presence of elastomeric damper in bearingless rotors, it is essential to 
characterise the non-linear behaviour of the elastomer and to develop a suitable model 
which can be easily integrated in the structural dynamic and aeroelastic analyses of 
the rotor blade. A review on the performance of elastomeric devices can be found 
in Ref. [22]. McGuire [23] conducted experiments at Lord Corporation to identify 
the non-linear behaviour of an elastomer. However, for subsequent analysis of the 
ground resonance problem, a linear mathematical model for the elastomer was used. 
Experiments carried out by Housmann [24, 25] indicated that most of the elastomeric 
materials show significant non-linear behaviour so far as the frequency and amplitude 
of motion are concerned. He modelled the elastomer in terms of strain-amplitude de- 
pendent complex moduli. In Ref. [24] he also proposed a thermo-mechanical model to 
include the temperature effects. Felker et al. [26] carried out experiments to determine 
the properties of an elastomeric lag damper used in the Bell Model 412 Helicopter. A 
non-linear model was proposed in which both the stiffness and damping characteris- 
tics were expressed as functions of displacement. Smith et al. [27, 28] extended the 
method of ADF [9] to include non-linearities and temperature effects. Recently, a 
model comprising of a series combination of a quartic spring and a linear Kelvin chain 
has been proposed by Gandhi et al. [29]. In this model, the values of the parame- 
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ters were determined by fitting the analytical curves with the experimental data. The 
non-linear hysteresis cycles at different equilibrium positions were examined. Later on, 
Gandhi and Chopra [30] proposed a modified elastomeric damper model to simulate 
a reduction in damping at very low dynamic amplitudes. With this improved elas- 
tomeric damper model, they showed a limit cycle oscillation for an autonomous system 
representing isolated lag dynamics of a blade. 

Panda and Mychalowycz [4] conducted experiments on a Comanche bearingless 
rotor with fiuidlastic (linear characteristics) and elastomeric (non-linear characteristics) 
dampers. However, in the aeroelastic analysis, the elastomeric damper w^as modelled 
by linear spring and damping elements placed between the flexbeam and the torque 
tube. Ormiston et al. [31] proposed a combination of a non-linear spring and a 
Kelvin chain with non-linear damping and spring elements. The damping element was 
represented by a linear combination of terms having fraction, linear, quadratic and 
cubic powers of velocity. In a recent publication, Kunz [32] discussed the merits and 
demerits of different modelling approaches for the elastomers, such as, ADF, complex 
modulus and VKS (Voigt-Kelvin-Solid). Based on the comparison of various models 
under harmonic excitation, he concluded that the non-linear VKS model would result 
in accurate estimates of energ\’ dissipation and was suitable for the prediction of forced 
response. However, this elastomer model was found to be amplitude-dependent. Using 
linear and non-linear VKS models for the elastomer, simulation studies of rotor start-up 
condition were compared with experimental results of AH-64A Apache main rotor. It 
was observed that neither the linear nor the non-linear model was capable of accurately 
predicting the response of the blade under a large impulsive loading. Of course, the 
non-linear model performed better than the linear one. The general conclusion of 
this study is that all elastomer models have certain limitations for applications to the 
analyses of the transient response, steady-state response and stability of rotor blades. 
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1.2.2 Structural Modelling of Rotor Blades 

Since helicopter rotor blades are long, slender beams undergoing moderate deforma- 
tions, a non-linear strain-displacement model is used to describe the coupling effects 
between axial, bending and torsional modes. Generally, the strains are assumed to be 
very small in comparison to unity. Such an assumption is consistent with the design 
requirement based on the fatigue life consideration which requires that the rotor blades 
must have an operating strain level well below the elastic limit of the blade material. 

The first structural model for a rotor blade undergoing flap-lag-torsional deforma- 
tions was developed by Houbolt and Brooks [33]. This model did not include the 
non-linear coupling effects between bending and torsion, which were shown to be im- 
portant by later researchers for the dynamic and aeroelastic analyses of helicopter rotor 
blades. The importance of these non-linear coupling effects has been discussed in a 
review article by Friedmann [34]. 

Depending on the method of formulation of non-linear beam theories, the structural 
model can be classified into two groups; (i) moderate deformation theories and (ii) large 
deformation theories. In the moderate deformation theories [35] - [40], higher order 
terms are elirpinated using an ordering scheme. In the large deformation theories, 
no such ordering scheme is employed. While the large deformation theories [41]- 
[44] are mathematically elegant, these are not easily amenable to aeroelastic studies. 
Therefore, moderate deflection theories are usually employed in the aeroelastic analysis 
of the helicopter rotor blades. 

Following the research and development efforts from late 70’s to early 90’s, the for- 
mulation of structural model of a rotor blade undergoing bending-torsion-axial modes 
has reached a level of maturity and is now more or less standardised. The procedure 
for the formulation of the structural model of a beam undergoing large (or moderate) 


8 



deformation consists of the following sequence of steps [45, 46]. 

i) The position vector of any point P in the undeformed configuration of the beam 
is defined with respect to a coordinate system rotating with the angular velocity of the 
rotor. 

ii) The undeformed base vectors at point P are evaluated. 

iii) The position vector of the point P is defined in the deformed configuration of 
the blade. 

iv) The deformed base vectors are then evaluated. 

v) Using the strain tensor in the curvilinear coordinate system, Green’s strain mea- 
sures are obtained. 

vi) The strains defined in the curvilinear system are then transformed into a local 
cartesian system. 

vii) In moderate deflection theories, the expressions pertaining to the axial strain at 
the elastic axis, the curv^ature of the deformed elastic axis and the twist of the beam are 
substituted with approximate non-linear expressions defined in terms of the physical 
displacements and slopes of the reference (elastic) axis. 

The moderate deformation beam theories developed by Hodges and Dowell [35] and 
Rosen and Friedmann [36] were validated by comparing the theoretical results with 
the experimental data obtained for the static deformation of a beam [40, 47]. These 
beam models were later used to formulate the inertial and aerodynamic operators for 
aeroelastic analysis of helicopter rotor blades [48] - [52]. In general, these models 
did not include the effects of cross-sectional warping and transverse shear, and were 
restricted to the treatment of isotropic blades. With the development of composite 
rotor blades, several researchers proposed structural models applicable for the analysis 
of composite rotor blades [53] - [59]. These composite blade models include cross- 
sectional warping and shear effects. A detailed review of structural models suitable for 
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composite rotor blades can be found in the literature [60] - [62]. 

While deriving the dynamic and aeroelastic equations of motion of a rotor blade, 
the inertia and aerodynamic operators generate a large number of higher order non- 
linear terms. A consistent ordering scheme is used to neglect all the higher order 
non-linear terms. The ordering scheme is based on assigning orders of magnitude to all 
nondimensional parameters of the problem, in terms of blade bending slope (which is 
assumed to be of the order of c « 0.10 to 0.15). A second order approximation assumes 
that terms of order are neglected in comparison to the terms of order 1. i.e., 

1 -f O K 1 

A few studies [38, 39] used a third order approximation where terms of order were 
neglected in comparison to terms of order 1. It must be borne in mind that the ordering 
scheme is based on experience and one should use certain degree of flexibility during 
its implementation. 

1.2.3 Stability and Response Analyses 

The problems in the area of stability and response analyses of a helicopter can be 
broadly divided into two categories, namely, (i) the problems related to the aeroelastic 
stability and response analysis of an isolated blade involving coupling between axial, 
bending and torsional modes; and (ii) the coupled rotor-fiiselage problems addressing 
the aeromechanical stability or flight mechanical stability or vibrations, depending on 
the range of frequency. In this thesis, the focus will be on the response analysis of an 
isolated bearingless rotor blade and aeromechanical stability (ground resonance) of a 
coupled rotor-fuselage system. In the following, a review of the pertinent literature is 
provided in two subsections each addressing a particular topic. 
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1.2. 3.1 Aeroelastic Stability and Response Analyses of an Isolated Rotor 
Blade 

An excellent review on aeroelastic stability and response of rotor blades can be found 
in [34]. Using moderate deflection beam theory, the aeroelastic stability analysis of 
an isotropic hingeless rotor blade was carried out by Hodges and Ormiston [48] and 
Shamie and Friedmann [49]. It was reported that a rotor blade configuration without 
precone is generally stable for all blade pitch angles; but a blade with precone can 
exhibit lead-lag instability for a variety of combinations of precone and pitch angles. 
The aeroelastic analysis of a hingeless rotor blade having a composite box beam in 
hover was presented by Hong and Chopra [63]. They observed that the coupling 
due to ply orientation has a strong influence on the stability of the blades. Later 
Panda and Chopra [64] extended this work to forward flight and demonstrated that 
aeroelastic stability could be improved by proper tailoring of the couplings generated 
by the composite lay-up. Recently Ttacy and Chopra [65] compared experimental and 
theoretical results in hover of a composite hingeless rotor blade model. 

The effect of tip sweep on the blade stability was first investigated by Celi and 
Friedmann [51]. They pointed out that tip sweep has a strong influence on the blade 
dynamics. Subsequently more studies have been carried out by different researchers 
[66, 67] to analyse the effects of swept tip on the rotor blade stability and response in 
both hover and forward flight. Yuan and FViedmann [46] made an analytical study on 
the aeroelasticity and structural optimization of composite hingeless blades with swept 
tip. It was observed that tip sweep has negligible influence on the fundamental flap 
and lag frequencies; but it can cause aeroelastic instability when there is a frequency 
coalescence between the first torsion and second flap modes. 

In a paper published in 1977, Bielawa [68] identified the complexities involved in 
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modelling a bearingless rotor blade. Three types of torque tube configurations were 
considered in this study. Two of the configurations exhibited performance and blade 
response characteristics which are similar to that of the conventional hingeless rotor 
blades. He concluded that from a practical design point of view “a snubbed torque 
tube configuration would probably be a minimum risk blade design” . A snubbed torque 
tube design is essentially the one shown in Fig. 1.3a, having an elastomeric snubber 
between the torque tube and the flexbeam. Sivaneri and Chopra [69] applied finite 
element method for the analysis of bearingless rotor blades. They have presented 
flutter stability of flap bending, lead-lag bending, and torsion in hover for a stiff-in- 
plane blade. They observed that the flap and torsion modes are stable while the lead-lag 
mode becomes unstable at high levels of thrust. This work was further extended by 
Hong and Chopra [70] for composite bearingless rotor blades. Dull and Chopra [71] 
investigated the effect of pitch-flap coupling for three pitch link configurations in a 
bearingless rotor blade. The hover stability predictions showed good correlation with 
the experimenatal data. They observed that the pitch link at leading edge provides 
higher lag damping compared to a pitch link at trailing edge. All these earlier studies 
[68] - [71] on bearingless rotor blades did not consider the non-linear elastomeric 

damper in their model. 

Recently, Gandhi and Chopra [30] developed a non-linear model for the elastomeric 
damper and integrated it with a blade model to examine the aeroelastic behaviour of 
a bearingless rotor blade. It was observed that due to the inclusion of the elastomeric 
damper, the fundamental lag frequency moved closer to the first rotor harmonic (1/rev.) 
resulting in an increase in the amplitude of periodic lag response. In a bearingless rotor, 
the pitch angle to the blade is provided by twisting the flexbeam by a control rod 
attached to the inboard end of the torque tube (Fig. 1.3a). Hence special attention 
must be paid to include the kinematic constraint at the inboard end of the torque 
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tube. The aeroelastic behaviour of the blade could also significantly depend on the 
initially deformed and twisted configuration of the flexbeam. Gandhi and Chopra 
[72] considered the effect of elastic twist of the flexbeam due to pitch link control to 
evaluate the aeroelastic stability characteristics of a bearingless rotor. It was observed 
that during the forward flight, for the same pitch setting of the blade, the leading edge 
pitch link configuration results in large positive torsional response of the blade (and 
hence an increase in the thrust) as compared to the trailing edge pitch link configuration 
which leads to a large negative torsional response (and hence a decrease in the thrust). 
It may be noted that in this study, the non-linear elastomer was not included. 

1.2. 3.2 Aeromechanical Stability of a Helicopter 

The aeromechanical stability of a helicopter is a complex phenomenon, involving cou- 
pling between the rotor and body degrees of freedom. The rotor lead-lag regressing 
mode usually couples with the body pitch and roll modes to cause instability. This 
phenomenon is denoted as ground resonance when the helicopter is on ground and air- 
resonance when the helicopter is in flight. The aeromechanical instability is severe in 
nature. In the case of articulated and hingeless rotor systems, the air/ground resonance 
is usually avoided by providing external mechanical dampers in the lead-lag mode of 
the rotor system, whereas in case of bearingless rotors, a specialised elastomer with 
high loss factor is used as a damper. Apart from providing an external lag damper, 
there are also other means of increasing the inherent damping available in the lag mode 
by judiciously providing certain aeroelastic couplings in the rotor blade. Another ap- 
proach is through the use of active controls based on fuselage and rotor feedback [73] - 
[77]. However, the safety and reliability of an active control scheme plays an important 
role which should be properly verified before its application to an actual helicopter. 

The analysis of aeromechanical stability chracteristics of a coupled' rotor/fuselage 
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system requires the development of a mathematical model representing the dynam- 
ics of the coupled system incorporating various aeroelastic couplings. Over the years, 
considerable progress has been made in modelling the system dynamics. Using a sim- 
ple mathematical model, Ormiston [78] examined the effects of pitch-lag and flap-lag 
couplings of a soft-in-plane hingeless rotor/fuselage system for both ground resonance 
and air resonance in hover. It was observed that aeroelastic couplings are generally 
not beneficial for ground resonance at zero-blade collective pitch. At higher collective 
pitch values, pitch-lag coupling is found to be more effective to stabilise the vehicle. In 
a later study, Ormiston [79] observed that the character of air-resonance was found to 
alter significantly at high rotor speeds for high lock number blade configurations. Fol- 
lowing the theoretical approach of Ref. [78], an experimental study was performed by 
Bousman [80] to investigate the influence of aeroelastic couplings and collective pitch 
setting on the aeromechanical stability of a hingeless rotor helicopter model. It was 
observed that the configuration with matched stiffness blade is basically less unstable 
than the configuration with non-matched stiffness. For the matched stiffness config- 
uration, addition of negative pitch-lag coupling (lag back-pitch up) was stabilising; 
however, for the non-matched stiffness case only slight improvement in the unstable 
region was observed. In general, the conclusions of Refs [78] and [80] indicate that the 
aeroelastic couplings ofier a significant potential for augmenting the inherent damping 
in the lag mode. 

Considering only the pitch-lag aeroelastic coupling, Zotto and Loewy [81] con- 
cluded that pitch-lag coupling could be effectively used to increase the damping in 
the lag mode under ground/air resonance conditions. Air resonance characteristics 
in forward flight with different dynamic inflow models have been studied by Nagab- 
hushanam and Goankar [82]. The study included the effects of trimming conditions, 
lag structural damping and aeroelastic couplings. Johnson [83], and Friedmann and 
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Venkatesan [84, 85] developed analytical models for the aeromechanical stability of a 
coupled rotor/fuselage system and compared their theoretical results with the experi- 
mental data presented in Ref. [80]. In Ref. [83], various cases considered for comparison 
between theory and experiment, were restricted to zero collective pitch setting of the 
rotor and in Refs. [84] and [85] the comparisons were limited to the cases without 
aeroelastic couplings. 

Venkatesan [86] analysed the effects of pitch-lag, pitch-flap and structural flap-lag 
couplings individually and in combination, on the stability of lag mode under ground 
resonance condition. He showed that the positive pitch-flap coupling increases the 
damping in resonance conditions, whereas the combination of the pitch-lag and flap-lag 
couplings do not have beneficial influence at resonance conditions. It was also observed 
that by a proper combination of flap-lag, negative pitch-lag and positive pitch-flap 
couplings, the lag regressing mode could be stabilised over the whole range of rotational 
speeds. However, this was done only at a single value of collective pitch. Recently, 
Gandhi and Hathway [87] extended the work of Venkatesan [86] by considering a 
combination of aeroelastic couplings that can stabilise ground resonance of a soft-in- 
plane hingeless rotor, over the entire range of rotational speeds, at various thrust levels 
and for various values of body frequencies. They considered the aeroelastic couplings 
as continuous variables and used optimization technique to reach a combination of 
couplings that completely stabilises the system under ground resonance. 

Gandhi and Chopra [88] analysed the ground/air resonance stability of a helicopter 
in hover with non-linear elastomeric lag damper. In their study, the blade was assumed 
to be rigid with a root hinge simulating both the articulated and hingeless rotor system. 

A non-linear elastomer was placed at the lag hinge and consequently it has the same 
dispalcement as that of the blade in lag mode. They observed that the non-linear 
elastomeric damper has a stabilising influence on both ground and air resonances. In 
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all the above studies [78] - [88], the blade model corresponds to a hingeless rotor 
blade. 

The dynamic analysis of a bearingless rotor blade is more complex than that of 
hingeless or articulated rotors owing to the presence of non-linear elastomeric damper 
and multiple load paths. One of the eeirliest analytical models for bearingless rotor 
was developed by Hodges [89, 90]. In the analysis, the main blade was assumed to be 
rigid. The flexbeam was treated as Euler-Bernoulli beam having elastic deformation 
in bending, torsion and axial directions. Based on finite element formulation of Ref. 
[69], Chopra and his associates examined the aeromechanical stability characteristics 
of a bearingless rotor in hover and forward flight conditions [91] - [93]. A study of 
the shaft-free aeromechanical stability characteristics of a bearingless rotor in hover 
[91] and in forward flight conditions [92] was carried out by Jang and Chopra. Tracy 
and Chopra [93] examined the aeromechanical stability of a bearingless rotor with 
composite flexbeam. This study indicated that a negative pitch-lag (lag back-pitch 
up) coupling has stabilizing effect on the regressive lag mode in hover and forward 
flight. In all these studies [89] - [93], the non-linear elastomeric damper was not 
included. 

Using the elastomer model proposed in Ref. [29], along with a rigid blade model 
having an elastomers at the root hinge, Gandhi and Chopra [88] examined its influence 
on the aermechanical behaviour in forward flight. It was observed that influence of 
damper dynamics on coupled regressive lag/fuselage motion could result in a diminished 
damping effect for some advance ratios. In another study, Gandhi and Chopra [72] 
also considered the effect of pitch link configuration on the shaft-free aeromechanical 
stability of a bearingless rotor without the elastomer. For the leading edge pitch link 
configuration, the regressive lag damping rapidly increases with advance ratios whereas 
for the trailing edge pitch link configuration, the regressive lag damping -decreases with 
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advance ratio. 

A bearingless rotor model was developed and tested by Weller [5] to examine the 
aeromechanical stability margins in hover. Rotor design parameters included funda- 
mental flap mode natural frequency, built-in cone and sweep angles of the blade, pitch 
link location at leading and trailing edge positions. It was observed that the stability 
characteristics at the higher collective pitch angles can be improved by attaching the 
pitch link to the leading edge. Recently, Panda and Mychalowycz [4] presented a 
correlation of wind tunnel tests and analytical studies for a scaled model of a RAH-66 
Comanche bearingless main rotor. However, in the analytical model the elastomer was 
treated as a linear element. Hover air-resonance test results were presented for both 
fluidlastic (linear charateristics) and elastomeric (non-linear characteristics) dampers. 
The experimental results indicated that the non-linearities in stiffness and damping 
seem to have a significant influence on the limit cycle oscillations. The authors at- 
tributed the reason for limit cycle to the presence of rotor-body frequency coalescence 
between the lag regressive, flap regressive and body roll modes. Near coalescence, the 
damping is low and hence a limit cycle oscillation occurs at the lag regressive frequency. 
On the other hand, the fluidlastic damper due to its linear nature did not exhibit any 
limit cycle oscillation. 

In Ref. [31], Ormiston et al. presented the capabilities of a Second Generation 
Helicopter Analysis System (2GCHAS). It is a finite element based comprehensive, mul- 
tidisciplinary computer programme. Using 2GCHAS, they studied the air-resonance 
stability characteristics of a 3-bladed hingeless rotor with rigid blades and a 5-bladed 
bearingless rotor with flexible blades. In the case of hingeless blade, a non-linear elas- 
tomer was placed at the root; in the bearingless blade, the elastomer was placed between 
the torque tube and the flexbeam. The elastomer was modelled as a non-linear Kelvin 
chain. The damping force was assumed to be proportional to a linear combination of 
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different powers of velocity (e.g., powers of 1/2, 1, 2 and 3). The phenomenon of limit 
cycle oscillation was demonstrated for both the configurations. It was concluded that 
“the effects of non-linear elastomeric dampers are significant and complex. Until the 
ideal bearingless rotor without auxiliary lead-lag damper is developed, the modelling 
and behaviour of elastomeric dampers warrant additional attention” . 

1.2.4 Non-linear Analysis 

Independent of the above studies related to helicopter blades, several researchers have 
made significant contributions to the study of non-linear dynamics of beams. The 
amount of literature available on the general non-linear problem is too vast to be 
considered in detail. In this section, the objective is to provide a review of the literature 
relevant to the present thesis. 

For small oscillations the response of an elastic body can be adequately described 
by linear, equations and linear boundary conditions. However, as the amplitude of 
oscillation increases, assumption of linearity does not hold good. In order to represent 
the dynamic motion, the non-linear effects have to be taken into account. The source of 
non-linearities may be geometric or material in nature. The geometric non-linearities 
may arise due to large deformation causing non-linear relationship between the strain 
and displacement. The material non-linearity occurs when the stress-strain relationship 
becomes non-linear in nature. 

FVom the experimental observation, it is well known that the dynamic stiffness and 
damping characteristics of a typical elastomer are highly non-linear with respect to the 
amplitude of deformation. In bearingless rotors, the presence of elastomer gives rise 
to material non-linearities. Additionally, moderate deformation theory produces geo- 
metric non-linear terms in the eqations of motion. In most of the non-linear dynamics 
problems, an exact solution is not available. Researchers generally follow either one 
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of the following schemes, namely, a) numerical techniques, b) approximate analytical 
techniques or, c) a combination of numerical-analytical methods. 

Numerical techniques deal with dynamic deformation using finite element or finite 
difference schemes in space and time. However, the application of the purely numer- 
ical techniques to vibration problems are costly in terms of computation time. On 
the other hand, analytical techniques are suitable for simple structures with simple 
boundary conditions. For complicated structures, majority of the researchers follow 
the numerical-analytic approaches. Using linear undamped natural modes, the govern- 
ing partial differential equations are converted into a system of coupled second order 
non-linear differential equations. These equations are solved by using perturbational 
technique such as method of averaging or method of multiple scales. This numerical- 
perturbational technique has been used for the determination of non-linear response of 
a continuous system with varying properties [94, 95]. 

Anderson [96] formulated the non-linear equations of motion of a rotating beam 
and obtained the natural frequencies. He concluded that the extensional and flapwise 
natural frequencies increase monotonically with rotational speed. However, his formu- 
lation contained certain errors which are pointed out in references [97, 98]. Non-linear 
flapping vibration of rotating blades had been studied by Venkatesan and Nagaraj 
[99] using the averaging technique. It was observed that the large amplitudes have a 
softening effect on the flapwise frequency. Experimental and theoretical analyses, of a 
composite blade under large deformation were carried out by Minguet and Dugundji 
[43, 44]. The analytical results w'ere obtained by using iterative finite difference solu- 
tion scheme. The natural frequencies and mode shapes obtained from the analytical 
model compared well with the experiment results. It was shown that equilibrium static 
deflection has a significant influence on the torsion and lead-lag modes and frequencies. 

Dowell [100] used the component mode analysis to examine the effect of material 
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non-linearity in the form of a non-linear spring-mass system attached to a simply 
supported beam. Pakdemirli and Naj’feh [101] extended the work of Dowell [100] 
by including the mid-plane stretching and damping. It was shown that both the mid- 
plane stretching and spring-mass system have a great effect on the frequency and force- 
response curves. Their relative effects depended on the system parameters and location 
of the spring-mass system. Birman [102] studied the free vibration of hinged beams 
on a non-linear elastic foundation. Szemplinska-Stupnicka [103] used generalised Ritz 
method to investigate the influence of non-linear boundary conditions on the frequency 
of beam vibration. Nayfeh and Nayfeh [104] obtained the non-linear normal modes and 
amplitude-frequency relation for a simply supported Euler-Bernoulli beam resting on 
an elastic foundation having quadratic and cubic non-linearities. Recently, non-linear 
normal mode shapes of a cantilever beam have been determined using the method of 
multiple time scale [105]. 

1.3 Objectives of the Present Study 

The review of the literature clearly indicates that the elastomer used in a bearingles 
rotor shows highly non-linear characteristics both in stiffness and damping with respect 
to the amplitude of deformation. Therefore, it is expected that due to the presence 
of the elastomer, the lag and flap frequencies should be dependent on the amplitude 
of motion which is a typical feature that distinguishes a non-linear system from a 
linear one. But so far no attempt has been made to study the effect of amplitude on 
frequencies of a bearingless blade. 

Several researchers have proposed different non-linear models with varying degree 
of complexity in order to capture the non-linear characteristics of the elastomer. How- 
ever, while performing the stability analysis, only the linearised equations of the blade 
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and elastomer about an equilibrium point are solved. In this solution procedure, the 
dependence of stiffness and damping characteristics of the elastomer on the amplitude 
of motion is completely ignored. 

All the studies on aeroelastic response of the rotor blade consider only the steady 
state response [32]. It was experimentally observed [106] that for a step input in pitch 
angle, the transient response of a hingeless rotor blade exhibits dynamic overshoot 
leading to overstressing of the blade. No information is available on the transient 
response characteristics of a bearingless rotor blade. Also no result is available in the 
open literature whether there is any qualitatively significant difference in the transient 
response of a bearingless rotor blade as compared to a hingeless blade. 

The objectives of the present study are to; 

i) formulate a simple non-linear model to represent the non-linear stiffness and 
damping characteristics of an elastomer; 

ii) study its influence on the isolated lag and flap dynamics of a rotor blade; 

iii) formulate a methodology for carrying out amplitude-dependent aeromechanical 
stability analysis under ground resonance condition considering only the elastomer 
non-linearity; 

iv) examine the transient response, to a step pitch input, of an isolated bearingless 
rotor blade undergoing coupled flapnlag-torsional deformation in hovering condition 
and compare the results with that of a hingeless blade; 

v) define a measure for the aeroelastic couplings in bearingless rotor blades, in terms 
of the steady state response (to a step pitch input) and study its dependence on the 
relevant geometric parameters. 

and 

vi) invesigate the phenomenon of limit cycle oscillation of an isolated bearingless 
rotor blade, with and without aerodynamic loading. 
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Figure 1.1: Articulated rotor hub and blade configuration. 
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Elastomer— 
Flexbeam — 




(b) Actual rotor hub (taken from Ref. [1]). 


Figure 1.3: Bearingless rotor hub and blade configuration'. 





Chapter 2 


ELASTOMERIC DAMPER 
MODEL 


2.1 Introduction 

In contrast to the conventional hydraulic dampers, elastomeric dampers contribute 
both damping and stiffness to the rotor blade. Due to the visco-elastic nature of the 
material, damping and stiffness properties of the elastomer are complex functions of 
the displacement amplitude, frequency and even temperature. A typical static load 
deflection curve for a non-linear visco-elastic material is shown in Fig. 2.1 by the curve 
OABC. Generally the static non-linear characteristics of the elastomer are described 
by the secant stiffnesses evaluated at different values of the deflection. For example, 
the secant stiffnesses at points A and B are defined by the slopes of the line OA and 
OB (Fig. 2.1), respectively. 

The load-deflection curve of the elastomer under dynamic loading about these two 
equilibrium points A and B can be represented by the corresponding visco-elastic 
hysteresis loops superimposed on the equilibrium curve OABC at the points A and 
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B, respectively. The dynamic stiffness values at A and B are defined by the slopes 
of the line pq and p'q' where the points p, q^ p' q' points correspond to 

the maximum deflection of the elastomer as shown in Fig. 2.1. On the other hand, 
linearised stiffness is defined using static load-deflection curve. The slopes of the lines 
uv and u'v' represent the linearised stiffness values in the vicinity of the equilibrium 
points A and B, respectively. 

In most cases, the dynamic tests on an elastomer is carried out without any static 
preloading i. e. about an equilibrium position with zero deformation. The material is 
subjected to a harmonic displacement excitation and the required force is measured. 
The stiffness modulus is defined as the in-phase force component that is required to 
produce unit amplitude of deformation. Similarly, the damping modulus is defined as 
the required force component (for unit amplitude of deformation) that is in quadrature 
with the deformation. The dynamic stiffness and damping values can also be identified 
from hysteresis cycles for different amplitudes of motion. The area of the hysteresis loop 
provides the information about the amount of energy dissipated per cycles at a par- 
ticular amplitude of deformation. Typical hysteresis cycles of a non-linear elastomeric 
material about zero equilibrium position are shown in Fig. 2.2. It can be observed 
from Fig. 2.2, that when the amplitude of motion is Xx, the corresponding dynamic 
stiffness is given by the slope of the line Oi6i. When amplitude of motion is increased 
to X 2 t the slope of the line 02^2 represents the dynamic stiffness. It is evident that 
these two stiffness values can be different depending on the distortion of the hysteresis 
loop with the amplitude of motion. However, the linearised stiffness about the zero 
equilibrium point is given by the slope of the static load deflection curve a 20 ,xOhxh 2 at 
0 i.e. slope of the line Im. It is important to recognise that the linearised stiffness 
can be different from the dynamic stiffness for a given amplitude. Also, the dynamic 
stiffness is a function of the amplitude of motion. Therefore in any linearised stability 
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analysis, the effect of change in the stiffness due to a change in the amplitude of motion 
is not at all taken into account. 

The energy dissipated by the viscoelastic material for a particular amplitude of 
motion is given by the area of the corresponding hysteresis loop. Figure 2.2 also 
indicates that the energy dissipation per cycle depends on the amplitude of motion. In 
case of a linearised analysis, a damping value has to be chosen based on an arbitrary 
value of amplitude of motion. 

The above discussion clearly indicates that the elastomer model should capture the 
non-linear behaviour of the material so far as the stiffness and damping are concerned. 
In addition, it should be possible to integrate the model easily with a structural dynamic 
model of a rotor blade for subsequent aeroelastic and aeromechanical analysis without 
neglecting the effect of the amplitude of motion on the elastomer characteristics. 

The experimental results of a single frequency bench test [23, 24, 26] of different 
elastomeric lag dampers are shown in Figs. 2.3-2.5. All the figures indicate that both 
the in-phase stiffness (denoted as K' or G') and the quadrature stiffness (denoted as 
K" or G") decrease as the amplitude of motion increases. The variation of loss tangent 
(tan (5) or loss factor ( 77 ) with amplitude is also shown in these figures. Loss tangent is 
generally defined as the ratio of G' and G". It can be noticed that there is no decrease 
in the quadrature stiiTness at very low amplitude of motion for all the experimental 
data shown in Figs. 2. 3-2.5. But in Ref. [30], Gandhi and Chopra assumed that 
there is a degradation of damping coefficient {G") at very low amplitude. This is not 
consistent with the experimental evidence shown in Figs. 2.3-2.5. However, as shown in 
Fig. 2.3b, the loss factor ( 77 ) decreases at low amplitudes. The reason for the decrease 
in loss factor at low amplitude of motion is mostly due to the sharp rise in the value 
of in-phase stiffness (G') at low amplitude. The computed loss factors for the other 
experimental data also exhibit similar behaviour as shown in Figs. 2.4b and 2 . 5 b. 
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In this chapter, an attempt is made to develop a model for the elastomer which 
would capture the amplitude dependent characteristics of stiffness and damping pa- 
rameters. A simple model is considered so that it can be easily integrated with a 
structural dynamic model of the rotor blade for subsequent dynamic and aeroelastic 
analysis. 

2.2 Formulation of the Model 

The non-linear characteristics of an elastomeric material can be idealised by a com- 
bination of linear and non-linear restoring and dissipative elements. Since the static 
stress-strain curve of an elastomeric material is non-linear, it is logical to assume the 
restoring element to be non-linear. The experimental results of a single frequency bench 
test [26] of an elastomeric lag damper show that both in-phase {G') and quadrature 
{G") stiffness decrease as the amplitude of motion increases; but neither of them dis- 
play any significant dependence on frequency within the range of interest (3.3 Hz to 6.4 
Hz). The experimental results taken from Ref. [26] are shown in Fig. 2.5, which are 
used in the present study. Figure 2.6 represents an idealised model consisting of a non- 
linear spring, a Coulomb damper and a hysteretic damper. The reason for choosing 
a Coulomb damper is that the damping force is very high at low amplitude. In addi- 
tion, Coulomb and hysteretic dampers provide damping forces which are independent 
of frequency. This type of idealisation differs from the model proposed by Gandhi et. 
al. [30], in the sense that their model has a viscous damping element. 

The non-linear force-deformation relation of the spring is assumed to be of the form: 

Fs = Kix- Kzx^ -i- RTgx® - Kyx’^ . (2.1) 

where Fs is the force exerted by the restoring element under a deformation x with ii'i, 
Kz, Kz and Ky as constants. The constitutive differential equation of the elastomer 
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model under harmonic loading (Fig. 2.6) can then be written as 


Kix - - Kyx"^ + Fsgn | x | +-x = Do sin u}t . 

CJ 


( 2 . 2 ) 


2.2.1 System Identification 

The parameters Ki, K 3 , K 5 , K 7 , F and h of the elastomer model are identified from ex- 
perimental data points, using error minimisation technique. The steady-state solution 
of Eq. 2.2 is assumed to be of the form 

X = X sin(a;t — /?) . (2.3) 

Substituting Eq. 2.3 in Eq. 2.2 and rewriting Do sinut as Do[sm(u}t— 0) coSy54-cos(u;t— 
P)sin 0 ], the following two equations are obtained by balancing the first harmonic. 

KiX - ^K^X^ + Ik^X^ - 1 ;K 7 X'^ = Do cos 0 (2.4) 

4 8 2 

4F 

and — + hX = DoSin0. (2.5) 

TT 

It may be noted that the first harmonic approximation of the function Fsgn | x | is 
given by [107] 

Fsgn I X I = F ^ cos (ut — 0) . 

Now, Eq. 2.3 can be rewritten as 

X = Xccosu}t + Xssincot (2.6) 

where Xc = —X sin 0 and Xs = X cos 0 . (2.7) 

Substituting for sin 0 and cos 0 in terms of AV and Xs in Eqs. 2.4 and 2.5, the in-phase 
stiffness (G') and the quadrature stiffness (G") are obtained, respectively, as 
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G' 


and G" 


Dc\Xs\ 

xi-^xl 

Ky - - \KrX^ 


Do\Xc\ 


' + ^1 


(2.8) 


(2.9) 


The experimentally obtained values for G' and G" for different amplitudes of motion 
are given in Fig. 2.5. Let these be denoted by G'{X) and G"{X), respectively, where A' 
is the amplitude of motion. The stiffness parameters JTi, Kz, amd Kj are obtained 
by minimising the error function: 

ND 

9, = EKWf (210) 

1=1 

where F, (A) = - ^K^X^ + ^ - ^ArA® - G" (A) , 

and ND denotes the number of experimental data points used. Similarly, the values 
of the damping parameters F and h, are obtained by minimising the error 

$2 

with Fj (A) 


ND 

E (Jf )]' 


( 2 . 11 ) 


2.3 Results and Discussion 

The experimental data is taken from Ref. [26]. These are given in FPS units which are 
converted to SI units for the present analysis. Number of data points ND is taken to be 
5. The frequency used for parameter identification is u; = 3.3 Hz. Gauss-Newton algo- 
rithm is used for the error minimisation problem. The system parameters so obtained 
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are given in Table 2.1. Figures 2.7 and 2.8 show the comparison of the experimental 
data with the theoretical results obtained using the idealised elastomer model. A fifth 
order approximation of the non-linear spring is also considered. From Fig. 2.7, it is 
evident that the seventh order approximation provides a better fit compared to the 
fifth order approximation. 

Since the hysteresis loop of a viscoelastic material provides information regarding 
the energy dissipation, an attempt is made to draw the hysteresis loop of the elastomer 
model about the zero equilibrium point. Figure 2.9 shows the hysteresis loops for 
different amplitudes of motion (X = 0.0005, 0.0010, 0.0015 and 0.0020 m). At the tip 
of the amplitude, the hysteresis cycles become almost vertical. This behaviour is due 
to the presence of the Coulomb damper in the model. 

2.4 Summary 

The most important aspects of the non-linear characteristics of an elastomers are sum- 
marised for convenience. The key difiierence between the linearised stiffness and dy- 
namic stiffness about an equilibrium point of the elastomer has been highlighted. Lin- 
earised stiffness is independent of amplitude of motion whereas dynamic stiffness is a 
function of the amplitude of motion. In addition, the damping of the elastomer is a 
function of the amplitude. 
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Deflection 


Fig. 2.1 Non-linear dynamic characteristics of a visco-elastic material. 



Fig. 2.2 Hysteresis cycles of a non-linear material about zero equilibrium point. 
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Fig. 2.3 Bifeenmental data taken from Ref. [24] 
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Fig. 2.4 Experimental data taken from Ref. [23] 
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Figure 2.6: Elastomer model. 




X X 10-^ 


Fig. 2.7 Variation of G' (in-phase stiffness) with amplitude, o experimental data; 
seventh order spring model; fifth order spring model. 



Fig. 2.8 Variation of G" (quadrature stiffness) with amplitude, o experimental data; 
— proposed model. 
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Force (N) 



Displacement (m) 

Fig. 2.9 Hysteresis loops of the elastomer model for different amplitudes of motion. 

X = 0.002m; X = 0.0015m; A” = 0.0010m; 

X = 0.0005m. 
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Table 2.1 System parameters of the elastomeric damper model. 


System 

Parameter 

Fifth Order 
Approximation 

Seventh Order 
Approximation 

2.543437 x 10” 
8.430654 X 10^^ 
1.084698 X lO^’' 

4.797347 x 10^ 
4.569120 X 10® 

2.673989 x 10® 
1.315287 X 10^2 
3.519586 X 10^^ 
3.176266 X 10^2 
4.797347 x 10^ 
4.569120 X 10® 

Ki (N/m) 
Kz (N/m^) 
Kz (N/m®) 
Kr (N/m^) 
F(N) 
h (N/m) 
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Chapter 3 


NON-LINEAR DYNAMICS OF 
UNCOUPLED FLAP AND LAG 
MOTION 


3.1 Introduction 

Dynamic analysis of rotating beams plays an important role in the design of vari- 
ous engineering systems such as turbo-machinery, robotic manipulators and helicopter 
blades. The natural frequencies and mode shapes of such structures have been a topic 
of constant interest and hence received considerable attention. In bearingless helicopter 
rotor, with the advancement in technology, the external hydraulic damper is replaced 
by an elastomeric damper. Experiments show that the stiffness and damping charac- 
teristics of an elastomer are highly non-linear in nature. Also the coupled equations 
of flap-lag-torsion and axial motions of a rotor blade are non-linear due to geometrical 
couplings. It may be noted that it is not known as to how much the non-linearities of 
the elastomer alter the dynamic characteristics of the rotor blade when the elastomer 
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is included as a sub-system. 

Due to the presence of the non-linear elastomer, the lag and flap frequencies should 
be dependent on the amplitude of motion which is a typical feature of a non-linear 
system. In this chapter, the influence of elastomer non-linearity on the undamped 
natural frequencies of a rotating blade has been studied. First a linear analysis is 
carried out for the uncoupled flap and lag modes and the results are validated using 
two different solution techniques. For the non-linear problem, the frequency-amplitude 
relationships in flap and lag modes are derived in a closed form, using numerical- 
perturbation technique. 


3.2 Bearingless Rotor Model 

A schematic diagram of the bearingless rotor is shown in Fig. 1.3a. Figures 3.1 a and 
3.1 b represent simplified models of the bearingless rotor under isolated flap and lag 
motions respectively. Since only the uncoupled flap and lag modes are considered, the 
rotor blade is idealised as a line element along the elastic (or reference) axis of the blade. 
The elastomeric damper is placed between the torque tube and the flexbeam to have 
adequate deformation and consequent lag damping. The torque tube is represented by 
a massless rigid link EC as shown in Fig. 3.1. The blade and flexbeam are assumed 
to be an Euler-Bernoulli beam with uniform cross-section and flexural rigidity. Since 
the damping does not play any significant role so far as the natural frequencies are 
concerned, elastomer is represented by a spring element. The effect of pitch link is not 
considered. Thus, for each mode, the bearingless rotor can be idealised as a rotating 
beam having a non-linear spring connected between the points B and C on the beam 
through the tome tube. A discussion on the free vibration analysis of the rotor blade 
system by assuming the spring to be first linear and then non-linear is presented in the 
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following.. 


3.3 Linear Analysis 

The dynamics of a rotating beam differs from that of a non-rotating one due to the 
addition of centrifugal stiffness. The differential equations of motion for a rotating 
uniform beam contain varaiable coefficients while those for a non-rotating uniform 
beam have constant coefficients. For the case of isolated flap dynamic problem, it is 
assumed that the elastomer is acting in the transverse direction only. Thus, there is a 
transverse constraint at the point B (Fig. 3.1 a) in the form of a linear spring whose 
deformation depends not only on the deflection at B, but also on the deflection and 
slope at the point C. On the other hand, in the case of lag dynamics, the constraint is 
in the in-plane direcion only (Fig. 3.1 b). 

In this section, the equations of motion and the associated boundary conditions 
for the uncoupled flap and lag motion of the rotating beam are presented separately. 
For comparative purposes, two different solution techniques, one based on power series 
expansion and the other based on the Rayleigh-Ritz principle are used to determine 
the natural frequencies and the mode shapes. 


3.3.1 Power Series 

Referring to Figs. 3.1 a and 3.1 b, the partial differential equations for small amplitude 
of vibration of the rotating uniform Euler-Bernoulli beam can be given by 


El 


d^Wi 

dx^ 


2 dx 


{L^-x^) 


dWi 


dx 


-m 




at* 


+ mQ^Wi 


for 0 < X < Li, 
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— m 


for Li < X < L 2 



—m 



+ mQ'^W2 




+ mCt^Wz 


for L 2 < X < L 


(3.1) 


where t denotes time, m denotes the mass per unit length, El the flexural rigidity, 
fl the angular velocity and Wi (i = 1,2,3) the deflection at difierent segments of the 
beam. In Eq. 3.1, the top terms in the right hand side apply to the out of plane motion 
whereas the bottom terms apply to the in-plane motion. (Note: For convenience, in 
this section, the same notations (W, W) have been used for both out of plane and 
in-plane motion.) 

Introducing the following dimensionless quantities 

= ^ = 1- ^ = ?- = fw = 1.2) , a = ^ (3.2) 

and assuming a harmonic solution of the form 


Wi (^, t) = ha (^) exp {vjjt) , 


Equations 3.1 reduce to 


a 


and a 


.d^Wi 

d^W2 

d^Wz 

d^^ 


1 d 

\(l 

2d^ 

<ii\ 

1 d 


2d^ 

s/ 

1 

H 

1 

1 d 


2d^ 

r 'nl 


-XWi 


-XW 2 

-XWz 


= 0 for o<e<ei, 

= 0 for < e < 6 (3.3) 
= 0 for 6 < ^ < 1. 


where A = 


for out of plane vibration, 
for in-plane vibration. 
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For the beam which is clamped at ^ = 0 and free at ^ = 1, the boundary conditions 


are 


„,.(0) = 0. M=0a„d^iM)=0, ^ = 0. 


(3.4) 


Consider the free body diagrams of the beam AD and the link EC, shown in Fig. 3.2, 
where V and M are respectively the force and the moment acting on EC and AD. 
Assuming the transverse/in-plane deflection to be positive along the Z/Y direction, 
the expression for V and M are given by 


V = KiA and M = 7(^-6) (3.5) 

where A =[Wc- (6 - 6) “ ^b] is the extension of the spring, 

with Wc = deflection of the beam at C, 

Wb = deflection of the beam at B, 

W'c — slope of the beam at C, 
and Ki = linear stiffness of the spring. 


In addition to the boundary conditions given by Eqs. 3.5, the following continuity 
conditions at the locations B and C have to be satisfied: 


W'l (fi) = w, (6) , 






^^2(6) - (6 -^i) 


de 

dW2{^2) 

d^ 


’ 

- (6) 




(3.6) 


w^2 (6) = VFg (^2) , 


m ( 6 ) dWz ( 6 ) 


If- It. -a 




di 


(6) - (& - «.) - W, ({,) 


W, (6) - (& - fi) 


d( 

d( 


- Wt (6) 


<PW,{h) 

de 




(3.7) 
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where K* = ^^^ 7 -- 

Equations (3.3) are linear ordinary differential equations with variable coefficients. 
Their solutions can be expressed as a power series in the independent variable ^ as 

00 

= ECtf*-' for 0<{<6 , 

A:=l 

00 

^2(0 = Ed*?*-* for 6<?<& (3.8) 

A;=l 

00 

and m {{) = E for 6 < f < 1 . 

A:=l 

Inserting Eqs. 3.8 into Eqs. 3.3 and equating the coeflBcients of the like powers of 
the following recurrence relations are obtained after some algebraic manipulations : 

2fl (k *1- 3) (At + 2) (Ac -f- 1) kC/^^^ — (A; “b 1) kCi^ + ^Ac (Ac — 1) — 2 ~ 0 

2a(A: + 3)(Ar + 2)(A: + l)A:E>;fc+ 4 - (A: + l)A:E»fc+ [it(fc-l)-2a;2] = 0(3.9) 

2a(fc + 3)(* + 2)(A: + l)A:£;jfc+4-(A: + l)A:E)t+ [lt(A:-l)-2a;2j^* = 0 

for Ac = 1,2, 3.... 

Similarly, inserting the expressions (3.8) into the boundary conditions (3.4), (3.6) and 
(3.7), one obtains 

00 00 

C, = 0 , C2 = 0 , '£Ek{k-l){k-2) = 0 , 5:E,(A:-1)(A:-2)(A:-3) = 0 ; 

jfc=i fc=i 

(3.10) 


oo 


ib=l k=l 

00 


k-l 


0 

i: c* (fc - 1) ({-^ - 5: ft (*: - 1) = 0 

k=l i=l 

00 00 

E C/b (* - 1) (* - 2) ff-’ - ft (* - 1) (t - 2) = 0 


(3.11) 


k=l 


k=l 
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k=:l 




Ar=l 


k=l 


k-l 


k-1 


-ZDt(k-l)(k-2)(i-3)(t‘' = 0 

k=l 


and 


A:=l fc=l 

E (fc - 1) fe*-" -Zb, (k-l) e*-= = 0 


A:=l 


Jt=l 


ZO,(k-l)(k-2) + *■• (f2 - 6) 


A:=l 




ifc-1 


A:=l 


E'D*e2‘‘’-te-«i)E^* 

*=1 fc=i 

oo 

E&(*-l)(*-2)f2*‘' = 0 


ZB>t(k-l)(k-2)(k-3) ft* - K‘ 


k=l 

00 
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E^*?2‘''-(6-6)E^* 


LA:=1 


A:=l 




A:~l 


A:=l 


ZE,(k-l){k-2)(k-3)ft* = 0. 

ifc=l 


(3.12) 


If the power series is truncated at the P-th term, then there are altogether 3P un- 
known coefficients. Prom the recurrence relations (3.9) and the boundary conditions 
(3.10), (3.11) and (3.12), there are 3P simultaneous linear homogeneous equations. 
For non-trivial solution, the determinant of the coefficient matrix must vanish. Thus, 
setting this determinant equal to zero one gets the frequency equation which is solved 
numerically for unknown uj. 

In order to obtain the mode shape, an additional normalising condition 


VF(1) = 1 


(3.13) 


is imposed. 
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3.3.2 Rayleigh-Ritz Method 

In this section, the equations of motion of the rotating blade under uncoupled flap and 
lag motion are separately derived using the Rayleigh-Ritz method. 


3. 3. 2.1 Equation of Motion in Flap Mode 

Figure 3.3 a shows the axial and out of plane deformations of a point A on the elastic 
axis due only to the flap bending motion of the rotating beam. The kinetic energy of 


the beam under is given by 




(3.14) 


Assuming that the blade is very stiff in the axial direction, the following inextensionality 


condition is imposed: 


“ = -5/ 
2 Jo 


1 fdw 


(3.15) 


The strain energy must include the effects of the bending moment, and spring defor- 
mation. The expression for strain energy due to bending is 




1 


The strain energy of the spring due to deformation A/ is 
V2 = 

where A; = w(L 2 ,t) - {L 2 - -w{Li,t) 

Hence the total strain energy of the system is 

U = Ui + U2. 


(3.16) 


(3.17) 

(3.18) 


(3.19) 


The equation of motiou is derived after introducing the dimensionless quantities defined 
by Eq. 3.2. Assuming a solution in the form 


n (r) 


(3.20) 
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motion. The above eigenvalue problem is solved numerically to determine the natural 
frequencies. 


3. 3. 2.2 Equation of Motion in Lag Mode 


Referring to Fig. 3.4, the kinetic energy of the beam under isolated lag motion is given 

by 

1 r r ^ { firi 1 ^1 

dx. (3.25) 


The expression for strain energy due to in-plane bending is 

d^v {x,ty^ 


Ui 


1 \i 
= - / El - 

2 Jo 


dx'^ 


dx 


(3.26) 


The strain energy of the spring due to deformation Aj in the lag direction is given by 


t/j = -ATiAf 


where A/ = 
Assuming a solution in the form 


s (i2, t) -(£,-£,) _ a 


dx 


(3.27) 

(3.28) 


= (3.29) 

»=i 

and following the same procedure as mentioned for the flap mode, the eigenvalue prob- 
lem for the blade in lag mode can be written as 


|K1{?} = A1M){,} (3.30) 


where [/f] = 


[M] = 



El 

& = 



Ki 

Ql = 

mQ^L 

and A = 
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3.4 Non-Linear Analysis 

The elastomer is represented as a non-linear spring. Following Pakdemirli and Nayfeh 
[101], a numerical-perturbation technique is adopted for the non-linear analysis of the 
rotor blade. Accordingly, the continuous system has to be first discretized by any one of 
the methods of weighted residuals [108] (such as. Collocation, Galerkin, Rayleigh-Ritz 
etc.) to obtain a system of coupled non-linear ordinary differential equations. Then, 
one can apply any one of the perturbation techniques. Though Galerkin discretisation 
procedure is most common, Ra,yleigh-Ritz technique is used in this study, since it gives 
a symmetric stiffness matrix for this problem whereas, the Galerkin procedure leads to 
a non-symmtric matrix. 

3.4.1 Non-Linear Equation of Motion 

First, the derivation of the non-linear equation of motion and close form expression for 
the frequency-amplitude relationship in flap mode is given in detail. For the lag mode, 
the seime procedure is followed. 

3.4.1. 1 Flap Motion 

When the elastomer is represented as a non-linear spring, the expression for the strain 
energy of the spring is modified as 

J/j = jATiAj - ijTsAj + gATsA' - iifjA; (3.31) 

Introducing the following dimensionless terms for the flap motion, 

X 

= — . {= - , r = at (3.32) 

with, 7 = j {<< 1), where f is the radius of gyration of the beam cross-section about 
the reference axis. Substituting for A/ in Eq. 3.31 from Eqs. 3.18 and using Eqs. 3.20 
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and 3.22, one can rewrite Eq. 3.31 as 


U2 = 




(t) 

,1 = 1 


4,2 


-KzL'^e 




Lt=l 




1 6 


U=1 


l 8 


1=1 


r. (r) 
(3.33) 


where e = 7 ® is the small parameter. 

Replacing the strain energy expression of the linear spring (Eq. 3.21) by the strain 
energy expression of the non-linear spring given in Eq. 3.33 and following the procedure 
outlined in Sec. 3.3.2. 1 , the temporal equation of motion in flap mode can be obtained 


as 


[M] {r} + [K] {r} -h {Fnl (r)} = 0 (3.34) 


where JF)vz (’’) contains the contribution of the non-linear terms in the spring force. 

In order to decouple the linear part of the above equation, the following linear trans- 
formation is introduced. 

r (t) = [P] z (r) (3.35) 


where [P] is formed by the normalised (with respect to mass matrix) eigenvector of the 
linear eigenvalue problem obtained after neglecting the non-linear terms in Eq. 3.34. 

Substituing Eq. 3.35 into Eq. 3.34 and premultiplying by [P]^, the temporal 
equation is obtained in the form 


n n n 

Zi (r) -I- ujfz, (r) - CQ;3 53 X) S ^imsrZmZsZr 

m=l 5=1 r=l 
n n n n n 

^ ^ ^ ^ ^ ^ ^ ^ ^ r* imsrtu^m^s^r^t^u 

m=l 5=1 r=l t=l ti=l 
n n n n n n n 

- X X X X X X X ^imsTtuvwZmZsZrZtZuZyZ^ = 0 (3.36) 

m=l 5=1 r=l t=l u=l v=l tiU=l 

i = l,2,3,...,n 


where u>i is ith linear natural frequency and 


tllsIWTRAL LrBIfAgl 

I- 1. T.. KAmmJi 


j ±1a£l>y 

"■ vTrXTwi 
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ri,„5r = tpi , 

r imsrtu = i’i , 

^tmsrtuvw ~ 

KzL KsL^ KrL^ 

mO,^ ’ mO^ ’ ’ 


(3.37) 


and 

= ^{pu + Up 2 i + - + (3.38) 

can be interpreted as the deformation of the spring in ith flap mode. In Eq. 3.38, 
Py represent the elements of the matrix [P]. The linear part of Eq. 3.36 is no longer 
coupled. Now, using the method of multiple scales and seeking an expansion of the 
solution of Eq. 3.36 for small but finite amplitude in the form 


(r ; c) = 2.0 (To, Ti , T2, Ta) + C2,i (To, Ti, T2, T3) + c"2i2 (To, Ti, T2, T3) + • • •• (3.39) 


where c is the small non-dimensional parameter. 
Different time scales are introduced as 


T„ = e„T , n = 0, 1,2, •••• (3.40) 

Here To is a fast scale associated with the changes ocurring with frequency close to Ui, 
while T„ for n > 1 are slow scales associated with changes that can only be noticed 
after several cycles. 

In terms of usual operator notations, we can write 

4- = T>o-}-£Di4-£^T»2 + c^T>3 + ---- (3.41) 

dr 

■^ = Dl + 2eDoDi -h (2D0D1 -f- D^) + { 2 D 0 D 3 + 2 D 1 D 2 ) + • • •• (3.42) 

where Di = , i = 0, 1, 2,* • •• 
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Substituting Eqs. 3.39 and 3.42 into Eq. 3.36 and equating the coefficients of like 
powers of e, one obtains 

: DflZio + = 0 (3.43) 

:Dlzii+J^Zii = -2DQDxZ^Q-\-OizTiiZ% (3.44) 

: DQZt2 <^i Zi2 = —2DoDiZii — D^Zio — 2DQD2Zio + SQ3T4iZ^QZ^i 

- ttsEeizfo (3-45) 

: DQZis + ujfzi3 = —2DoDiZi2 — D^zn — 2DoD2Zii — 2DoD3ZiQ — 2Z?i 1 ) 22,0 

+3a:3r4i2?o2^i2 + 3Q:3r4,2i0-Z,^l “ 5Q;5r6.2,V»i + aTEs.^Jo (3-46) 

where 

Ui =-{*{)' , r6. = (,t/)‘ , r„ = (^/)® , (3.47) 

Assuming the solution of Eq. 3.43 as 

Zio = .4j(Ti,T2,T3)e"'’'" + ii(r,,r2,r3)e-"‘^« . (3.48) 

Substitution of the first order solution given by Eq. 3.48 into Eq. 3.44, yields 

Dlza + u;^zn = -2iuJiDi (Ai) + a3r4i + 3AjAre^'^°] + c.c. (3.49) 

where c.c. indicates the complex conjugate. 

In order that zulzio be bounded for all To , the secular terms must vanish. Hence 

— 2iu)iDi (Aj) + 3Q;3r4jA|.(4j = 0 (3.50) 

and the solution for zn becomes 

Zii = -^a 3 r 4 iAfe*^'^“ + c.c. (3.51) 

ouf 

In order to uniquely define the amplitude A, of the fundamental frequency of oscillation, 
the homogeneous solution of Eq. 3.49 is not included. 
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To solve Eq. 3.50, one can assume 


Ai = -a,e‘>‘ 


(3.52) 


where Oj and /?, are real functions of Ti, T2, T3. Substituting equation (3.52) into 
equation (3.50) and separating the real and imaginary parts, one obtains 

do, d0i 3 ^ 2 

Hence a, = , ft = -j^Qar.iaJr, + ;»; (Tj.T,) . (3.53) 

O COi 

Substitution of Eqs. 3.48 and 3.51 into Eq. 3.45, one gets 
Dlzi2 + Ujfzi2 


- 2iuiD2 (Ai) 


32a;? ^ 64a;2 


8 a;? 3 '' ' '' 8 a;2 " 4 a;? 

-OsFei + c.c. 


(3.54) 


Again the absence of any secular term in the solution of Eq. 3.54 requires 
— - 2iuiD2 (Ai) - ^alrlAfA^ - lOosEeiAfi? = 0. 


(3.55) 


-- 8^2 

Separating the real and imaginary parts, we have, after some algebraic manipulations, 

n 9 - 2r>2 , , . „ 3 _2r2 ^5 3 


dT 2 ^ ’ 128 a;? 


+ <^iai 256 a;2 ^^OisTeia, - 0. 


15 5 

Hence a = a, (T,) ./j: = - jgg^alrJiafr, + ^asHMafT, + /3" (T,) . (3.56) 

The solution of Eq. 3.54 is given by 

- __?Z_rv2r? -4- ^ n?r? i** J Ow.To , ^ ^,2p2 46 tSw.To 

" 256 a;4 A ^ ^4 "si 4, A Ae + ^^2 ^4 "s^ 4i>ii e 


.^a5re,4e'=“'^" + g^a5r6i.4Mie‘^''" + c.c. 


(3.57) 
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Substitution of Eqs. 3.48, 3.51 and 3.57 into Eq. 3.46, yields 


DqZiz + Ziz — —2 DqDi 


._^Z_ 2p2 243 t3a>.ro , _i_^2p2 ^4 J i3a;.7b 

256 "^3^ 4, -a. A e +32 ^4 0 = 3 ! 4^A Ae 




24a;f 


-i.—a5reiAfA,e'^^^^ + c.c. 


-2D0D2 + c.c.j - Dl + 

-2DoDz [A,e^'^° + c.c.) - 2 D 1 D 2 + c.c.) 


+3a3r4. + C.C.) ^-^a3r4,A'e’^-^° + c.c. j 

-5a5r6. + 6Afi2 + c.c.) 

+3a3r4, ^^2gt2a;.ro 2AA + C.C.) 2,2 + OtFs, ^^, 6“^’^® + c.c.)^. (3.58 


The condition for eliminating the secular term can be given, after some algebraic ma- 


nipulations, as 


. aa, 80'^ 5 _ _ 7 111 3T.3 7 35 ^ 7 ^ 


(3.59) 


Again, separating the real and imaginary parts of Eq. 3.59, one obtains 
9at „ _ . . 


= 0 or, a, = constant 


^80" 5 ^ ^ 6 111 , , 6 

&Tz ~ 64 x 128a;,5“"^'°‘ 

// 5 111 

CA , .3*^3r'4tQ:5r6iaj T 3 — TOO. .5^3^4i^t 


»T> ~ CA 3“3^ 4t'^5^ 6t«'i s'^3'^ 4t“i Too 

8Tz 640;; 64 x 128a;f 128o;, 


arEsiaf 


or, 0i = 


64 X 128 a;? 


c^TA^a^Tz 


128 a;,- 


Ol7^Si<liTz A 00 = 0 


(3.60) 


Prom Eqs. 3.40, 3.53, 3.56 and 3.60 one can write the frequency-amplitude relation as 
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+ e- 


ra3r4,-Qr5r6taf — 


111 


35 


«|ri.af - e^^a-rFsia^ (3.61) 


64^3--^- — - 64 ^ ' 128a;, 


where uj^i is the e-th non-linear natural frequency. Redefining the amplitude as 




Qi 

V~e 


(3.62) 


the final frequency-amplitude relation can be written as 


3 ■r\ A 9t^9 j ^ ^ T-\ j ^ 

® a3r4ia5re.ii‘ - - —a-rV^iAt. (3.63) 


64 0; 


64 X 128a;f 


128 o^i 


3.4.1. 2 Lag Motion 


The elastomer is represented by a non-linear spring. It is assumed to be effective in 
the in-plane direction only. The assumed mode solution in lag mode is given by 

»=i 

where is the i-th rotaing mode shape in lag motion. The expression for can be 
written as 

+ ^2P2i + ... + ^nPn, (3.64) 


with pij as the elements of the matrix formed by the eigenvetor of the linear eigenvalue 
problem. 

The expression for the strain energy of the spring in nondimensional form can be 
written as 




w 


U=l 
n 


- -K^L^e^ 


U=1 




,1=1 




n8 


Y1 


Li=l 


(3.65) 


where i/)- can be interpreted as the deformation of the spring in the lag mode. 


56 



Following the same procedure as outlined in the earlier section on flap vibration, 
the non-linear temporal equation for lag vibration can be obtained which will have the 
same form as given in Eq. 3.36. Using multiple-time-scales, the final expression for 
amplitude-frequency relation for the isolated lag mode of the rotating beam is obtained 
as 




3 p : 2 

Ui - 0:3r4i4 

8uji 




256 






5 -'6 

r 0^3 f 4i ^5 .^4, 


1X1 o_,o ^ 6 


64 w? 64 X 128 a.-. 

where Wj is ith linear natural frequency in the lag mode and 


OjTsiA" (3.66) 


Fimsr = . 

Fim^rtu = i>\‘4>'rni>Wr'^tipu - 

KzL K7L^ 

mQ2 ’ - ^^2 ’ “7 - ^^2 ■ 


( 3 . 67 ) 


3.5 Results and Discussion 

This section describes the numerical results obtained for the linear and non-linear 
dynamics of the blade in uncoupled flap and lag motions. First, the results pertaining 
to the rotating beam in the flap mode are presented. 

3.5.1 Isolated Flap Motion 

The results for the out of plane vibration of the rotating beam are given below in three 
parts. 
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3. 5. 1.1 Rotating Beam Without Elastomer 


The beam is assumed to be uniform and the data used in the computations are given 
in Table 3.1. These values are taken from Refs. [109, 110] for correlation purpose. The 
natural frequencies and mode shapes of a rotating beam are obtained by using both 
the power series and the Rayleigh-Ritz method. The first three natural frequencies 
in the flap mode are calculated. The roots of the frequency equation obtained by the 
power series method are computed by an iterative search procedure. For the Rayleigh- 
Ritz analysis, the following comparison functions [111] satisfying all the boundary 
conditions are used : 

■A. (0 = e 

where i refers to the mode number. 

The results so obtained are given in Table 3.2. It is observed that the results 
obtained by both the methods are in excellent agreement with those presented in Refs. 
[109, 110]. It may be noted that in Refs. [109, 110], the natural frequencies have been 
obtained by finite element analysis of the blade. 

3.5. 1.2 Rotating Beam with Linear Elastomer 

The elastomer is idealised as a linear spring. The dynamic characteristics of a rotating 
beam are analysed for different values of elastomer (=^) and torque tube ^2 (=^) 
locations (see Fig. 3.1 a). Calculations are also performed for four values of the dimen- 
sionless spring constant K* =1000, 2000, 3000 and 4000. To validate the results, the 
frequenceies are evaluated by using both the power series method and Rayleigh-Ritz 
technique. The power series converges for P — 50, where P represents the number of 
terms in each series. Since the rotating beam is divided into three parts and each part 
is represented by a 50 term power series, there are in total 150 unknown coefficients. 


(i -f- 2) {i -t- 3) ^ i 3) 2 


-e- 


r + 


* (* + 1)^3 


(3.68) 
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While using the Rayleigh-Ritz technique, the natural frequencies converged when six- 
teen comparison functions {i =1 to 16 in equation (3.68)) were used. Table 3.3 presents 
the first three natural frequencies of the rotating beam, evaluated by both techniques, 
for different values of the parameters [K*, la)- In most of the cases, the natural fre- 
quencies obtained by the two methods are in excellent agreement with each other. 
However, there is a difference of about 3% to 6% in the third natural frequency for 
la = 0.3. For a given configuration, say |i = 0.1 and |a = 0.25, when the stiffness (K*) 
is wied from 1000 to 4000, the percentage variation in natural frequency in the first 
mode is 0.6%; in the second mode the variation is 0.02% and in the third mode the 
\’ariation is 2.17%. Similar behaviour is also observed for other values of |a. 

The mode shapes corresponding to the first three natural frequencies are shown 
in Figs. 3.4-3. 6, for |^i = 0.10, la = 0.20, 0.25 and 0.30 and jFf*=2000. The mode 
shapes of the rotating beam without the spring are also included in these figures. It 
is observed that the presence of the linear spring affects the mode shape in the region 
between |i = .10 to la = .40. The influence of the spring is more pronounced for the 
third mode (Fig. 3.6). 

A detailed analysis is carried out to study the influence of different geometric con- 
figurations of the torque-tube on the natural frequencies. The value of |i was varied 
from 0.11 to 0.25 in steps of 0.02 and that of |a is varied from 0.25 to 1.0 in steps of 
0.05. The results are generated for academic interest only since no practical blade will 
have la > 0.3. The value of the spring stiffness was chosen as K* = 2000. The natural 
frequencies are evaluated by the power series method. The variation of the first three 
natural frequencies are shown in Figs. 3. 7-3.9. From Fig. 3.7, it can be seen that for 
all values of |i, the first natural frequnecy initially increases and then decreases with 
increasing |a. The maximum value of the frequency occurs around |2=0.40. The sec- 
ond natural frequency also exhibits a similar trend (Fig. 3.8) with the exception that 


59 



the maximum value takes place around ^2 =0.70. The third natural frequency shows 
more complicated nature of variation with increasing ^2 (Fig. 3.9). It may be noted 
that the point of attachment of the torque tube (i.e. ^ 2 )? rather than the location (^ 1 ) 
and stiffness {K*) of the spring, predominantly governs the natural frequencies. 

3. 5. 1.3 Non-Linear Analysis 

The force-deformation characteristics of the non-linear spring are represented by two 
models: one with a fifth order approximation and the other with a seventh order 
approximation as detailed in Chapter 2. The fifth and seventh order approximations 
are given as 

Fs — Kix — and Fs = Kix — K^x^ -1- K^x^ — K 7 X^, respectively. 

The values of the constants are given in Table 2.1. For the fifth order approximation, 
corresponding to the value of Ki, given in Table 2.1, is K*=3480A. Similarly for the 
seventh order appproximation, JT*=3658.8. It should be noted that both these values 
are close to 4000 for which the numerical results with the linear analysis have already 
been presented. The non-linear dynamic analysis is carried out for a specific rotor 
blade configuration having = 0.1 and ^2 = 0.25. These are representative values 
corresponding to a realistic bearingless rotor. 

The procedure of evaluating the frequency-amplitude (Eq. 3.63) relation is as fol- 
lows. Using the linear stiffness term from the nonlinear force-deformation relation of 
the spring, a linear eigenvalue problem is solved for the natural frequencies, the trans- 
formation matrix [P] (Eq. 3.35). [P], the deformation of the spring in each mode ( 
■ipl ) is obtained from Eq. 3.38. Knowing a,’s from Eq. 3.37 and using Eq. 3.47, the 
coefficients of the frequency-amplitude Ekj. 3.63 are obtained. The final expressions of 
the frequency-amplitude relationships for the first three modes corresponding to the 
two spring models, are given below: 
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5th order spring model : 


ujm = 1.138 - 20.1568^2 - 1.487626454 x 10l4'‘ + 3.188617261 x 

wni = 3.408 - 3.947019815 x - 1.904996008 x 10“^i^ + 4.780973363i'‘ 

ujm = 7.924 - 8.797674851 x loM^ - 4.070122676 x lO^i'^ + 2.426085485 x 10^°i'‘ 

7th order spring model : 

u)Ni = 1.138 - 2.793^2 - 2.8553 x lO^i^ + 8.657906 x lOM'* + 1.41639283 x 10®i® 

- 4.32084 X lO^i® - 9.50547154 x lO^i® 

UJM 2 = 3.408 - 6.7618 x lO^^i^ - 5.5908 x lO'^i^ + 17.850li^ + 2.36143 x lO'M® 

- 6.8414 X 10"^i® - 1.66869 x lO^i® 

ujNz = 7.930 - 1.226218 x lOM^ - 7.900145i^ + 6.650156i'‘ + 6.855182 x 10^''i® 

- 7.532292 x lO^^i® - 1.277168 x 10^®i® 

The variations of the first three natural frequencies with the amplitude of oscillation 
at the free end are shown in Figs. 3.10-3.12. From Fig. 3.10, it can be seen that with the 
seventh order approximation, the first natural frequency decreases continuously with 
the amplitude. Where as, with the fifth order approximation, the frequency initially 
decerases upto an amplitude of 0.018 and then shows an increase with a further increase 
in the amplitude. The reason for this behaviour can be attributed to the stiffening 
nature of the fifth order approximation of the spring at high amplitudes, as shown in 
Figure 2.7. A similar behaviour is depicted by the third natural frequency as can be 
seen from Fig. 3.12. The second natural frequency is least affected by the amplitude 
upto a value of 0.1 as can be seen in Fig. 3.11. 

Figure 3.13 shows the nature of variation of the non-linear first natural frequency 
with amplitude as one uses different orders of approximation for the spring. With 
the inclusion of higher order non-linear term (one at a time), the frequency-amplitude 
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relationship shows alternate softening and stiffening nature. Therefore, the range of 
validity of the frequency-amplitude relationship (for a particular order) is limited to 
the point beyond which there is appreciable variation. 

3.5.2 Isolated Lag Motion 

Numerical results corresponding to uncoupled lag motion are obtained using the method 
outlined in Secs. 3.3 and 3.4 and presented in three parts, each one concentrating on 
a particular problem. The data used are given in Table 3.1. 

3.5.2. 1 Rotating Beam without Elastomer 

The natural frequencies and mode shapes of a rotating beam in the lag mode are 
obtained by using both the power series and the Rayleigh-Ritz method. Since the 
third natural frequency is usually very high for the lag mode, the analysis is restricted 
to the first two frequencies only. For the Rayleigh-Ritz analysis, the same comparison 
functions given in Eq. 3.68 satisfying all the boundary conditions are used. The results 
are given in Table 3.4. It is observed that the results obtained by both the methods 
are in excellent agreement with those presented in Refs. [109, 110]. 

3.5.2.2 Rotating Beam with Elastomer 

The elastomer is idealised as a linear spring. Calculations are performed for three 
values of the dimensionless spring constant K* = 1000, 2000 and 3000. The dynamic 
characteristics of a rotating beam are analysed for different values of (=^) and ^2 
(=^) (Fig. 3.2 b). For the in-plane vibration, the power series converges for P = 30, 
where P represents the number of terms in each series. Since the rotating beam is 
divided into three parts and each part is represented by a 30 term power series, there 
are in total 90 unknown coefficients. While using the Rayleigh-Ritz technique, the 
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natural frequencies converged when sixteen comparison functions (i =1 to 16 in Eq. 
3.68) were used. The first two natural frequencies of the rotating beam, evaluated by 
both techniques, for different values of the pmameters {K*, ^ 2 ) aie included in Table3.5. 
In most cases, the natural frequencies obtained by the two methods are in excellent 
agreement with each other. However, there is a difference of about 2.4% in the first 
natural frequency for ^2 = 0.3. The difference is less (about 0.17%) for the second 
natural frequency. The results indicate that the point of attachment of the torque 
tube (i.e. ^ 2 ) rather than the stiffness {K*) of the spring, predominantly governs the 
natural frequencies. 

The mode shapes corresponding to the first two natural frequencies are shown in 
Figs. 3.14 and 3.15, for the values ^ = 0.10, 6 = 0.25 and F!:*=2000. The mode 
shapes of the rotating beam without the spring are also included in these figures. It 
is observ'ed that the influence of the spring stiffness (with K* — 2000) is not much 
pronounced. 

As in the case of the out of plane vibration, a detailed analysis is carried out 
to study the influence of different geometric configurations of the torque-tube on the 
natural frequencies. The value of the spring stiffness was chosen as K* = 2000. The 
natural frequencies are evaluated by the power series method. The variation of the first 
and second natural frequencies are shown in Figs. 3.16 and 3.17. Prom Fig. 3.16, it 
can be seen that for all values of ^ 1 , the first natural frequnecy initially increases and 
then decreases with incretising ^ 2 - The maximum value of the frequency occurs around 
^2=0.40. The second natural frequency also exhibits a similar trend (Fig. 3.17) with 
the exception that the maximum value takes place around ^2 =0.70. 



3. 5. 2. 3 Non-Linear Analysis 


The procedure of evaluating the frequency-amplitude relation of the flap vibration is 
mentioned in Sec. 3. 5. 1.3. Following the same procedure, the variation of the first two 
lag frequencies with amplitude of oscillation are obtained and are shown in Figs. 3.18 
and 3.19. These results indicate that the frequency-amplitude behaviour in lag mode 
is similar to that obtained for flap mode. 


3.6 Summary 

Both linear and non-linear free vibrations of a rotating blade under isolated flap and lag 
bending have been studied. For the linear problem, two different solution techniques, 
one based on the power series expansion and the other bsed on Rayleigh-Ritz principle, 
were used. Natural frequencies obtained from these two methods are found to be in 
excellent agreement for both types of beam motion. 

In the non-linear analysis, a numerical-perturbation technique is applied to de- 
termine the frequency-amplitude relationship for isolated flap and lag motions of the 
rotating beam. It is concluded that up to a fairly high-value of amplitude, a seventh 
order expansion is sufficient to correctly predict the frequency- amplitude relationship 
for both cases. If the amplitude level is lower, one can use an expansion up to fifth 
order and still get accurate results while saving on the computational effort. The al- 
ternating signs in the polynomial expressing the restoring force prevents the non-linear 
characteristic to be monotonically softening or hardening. 
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E 



Fig. 3.2 Free-body diagrams of the bearingless rotor blade and torque tube. 
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z 



(a) Flap 



(b) Lag 


Fig. 3.3 Elastic deformation of the blade. 
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Beam deflection 


Fig. 3.4 First mode shape of the rotating beam in flap mode for various values of ^ 2 - 


o, no spring. ^2=0.20; ^2=0.25; ^2 =-30; 

^1=0.10, K*=2000. 



Fig. 3.5 Second mode shape of the rotating beam in flap mode for various values of 


^ 2 - o, no spring. 6=0-20; 

6=0.10, K*=2000. 


6=0.25; 6 =0.30; 





Fig. 3.6 Third mode shape of the rotating beam in flap mode for various values of ^ 2 - 

o, no spring. ^2=0.20; ^2=0.25; ^2 =0.30; 

6=0.10, K*=2000. 



Fig. 3.7 Influence of different geometric configurations (6 and 6) on the first natural 
frequency in flap mode. iT*=2000. 
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Fig. 3.8 Influence of diiTerent geometric configurations (^i and ^ 2 ) on the second 
natural frequency in flap mode. A"* =2000. 



Fig. 3.9 Influence of different geometric configurations (^1 and ^ 2 ) on the third naturaJ 
frequency in flap mode. A* =2000. 
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Fig. 3.10 Variation of the first non-linear frequency with amplitude of oscillation in 
flap mode. 

fifth order approximation; seventh order approximation. 



Fig. 3.11 Variation of the second non-linear frequency with amplitude of oscillation 


in flap mode. 

fifth order approximation; 
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seventh order approximation. 






<^N3 

Fig. 3.12 Variation of the third non-linear frequency with amplitude of oscillation in 
flap mode. 

fifth order approximation; — seventh order approximation. 



Fig. 3.13 Influence of different orders of approximation on the first non-linear natural 
frequency in flap mode. 

third order approximation; fifth order approximation; seventh 

order approximation. 
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Beam defli 


Fig. 3.14 First mode shape of the rotating beam in lag mode 

o without spring. with spring. 

6 = 0.10, £2 = 0.25 and K* = 2000. 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


Fig. 3.15 Secone mode shape of the rotating beam in lag mode 

o without spring. with spring. 

6 = 0.10, 6 = 0.25 and K* = 2000. 







Fig. 3.16 Influence of diff'erent geometric configurations (^i and ^ 2 ) on the first natural 
frequency in lag mode. /sr*=2000. 



Fig. 3.17 Influence of different geometric configurations (^1 and ^ 2 ) on the second 
natural frequency in lag mode. K*=2000. 
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Fig. 3.18 Variation of the first non-linear frequency with amplitude of oscillation in 
lag mode. fifth order approximation; seventh order approximation. 



Fig. 3.19 Variation of the second non-linear frequency with amplitude of oscillation 
in lag mode. fifth order approximation; — seventh order approximation. 
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Table 3.1 Blade data. 


m 

L 

El 

mn^L* 

(Flap) 

El 

(Lag) 

r 

n 

9.7Kglm 

6.6 m 

0.0106 

0.0301 

.5 m 

32.8 rad f sec. 


Table 3.2 Natural frequencies of the rotating blade in flap mode without elastomer. 


Mode 


Rayleigh- 

Ritz 

Ref. [109] 

Ref. [no] 

1st 



1.125 

1.1247 

2nd 



« 


3rd 

7.6218 


- 

7.6376 


Table 3.3. Natural frequencies of the rotating blade in flap mode with a linear elas- 
tomer. 




6 = 0.1, 6 = 0.15 

= 0.1, 6 = 0.2 

6 = 0.1, 6 = 0.25 

6 = 0.1, 6 = 0.3 

pm 


Power 

Series 

Rayleigh- 

Ritz 

Power 

Series 

Rayleigh- 

Ritz 

Power 

Series 

Rayleigh- 

Ritz 

Power 

Series 

Rayleigh- 

Ritz 



1.125 

1.125 

1.127 

1.127 

1.130 

1.131 

1.133 

1.134 



3.407 

3.407 

3.409 

3.409 

3.407 

3.407 

3.413 

3.416 



7.617 

7.622 

7.628 

7.635 

7.716 

7.763 

7.960 

8.152 



1.125 

1.125 

1.128 

1.129 

1.132 

1.134 

1.134 

1.137 



3.408 

3.408 

3.409 

3.409 

3.407 

3.407 

3.414 

3.420 



7.617 

7.635 

7.634 

7.645 

7.754 

7.847 

8.039 

8.349 

1st 


1.125 

1.125 

1.130 

1.131 

1.133 

1.136 

1.135 


2nd 

3000 

3.409 

3.408 

3.410 

3.411 

3.407 

3.407 

3.415 


3rd 


7.618 

7.622 

7.639 

7.656 

7.774 

7.903 

8.075 


1st 


1.125 

1.125 

1.130 

1.133 

1.134 

1.138 

1.136 

1.139 

2nd 

4000 

3.409 

3.408 

3.410 

3.412 

3.407 

3.408 

3.415 

3.424 

3rd 


7.618 

7.622 

7.642 

7.664 

7.786 

7.937 

8.095 

8.509 
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Table 3.4 Natural frequencies of the rotating blade in lag mode without elastomer. 


Mode 

Power 

Series 

Rayleigh- 

Ritz 

Ref. (109] 

Ref. [110] 

1st 

2nd 

0.7317 

4.4825 

0.7317 

4.4825 

0.732 

0.7326 

4.4563 


Table 3.5. Natural frequencies of the rotating blade in lag mode with a linear elastomer. 



■■1 

6 = 0.1, 6 = 0.15 

6 = 0.1, 6 = 0.2 

6 = 0.1, 6 = 0.25 

6 = 0.1, ^2 = 0.3 

Mode 


Power 

Series 

Rayleigh- 

Ritz 

Power 

Series 

Rayleigh- 

Ritz 

. 

Power 

Series 

Rayleigh- 

Ritz 

Power 

Series 

Rayleigh- 

Ritz 

1st 

2nd 


0.733 

4.484 

0.733 

4.484 

0.744 

4.487 

0.745 

4.488 


0.768 

4.484 


0.792 

4.493 

1st 

2nd 


0.734 

4.485 

0.734 

4.485 

iiil 

0.755 

4.492 

0.774 

4.484 

0.786 

4.485 


^ 0.808 
4.497 

IBM 


0.735 

4.486 

0.735 

4.486 

0.756 

4.492 

0.764 

4.496 


0.796 

4.485 

0.796 

4.491 

0.815 

4.499 
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Chapter 4 


AMPLITUDE DEPENDENT 
STABILITY OF AN IDEALISED 
BEARINGLESS ROTOR UNDER 
GROUND RESONANCE 

4.1 Introduction 

The aeromechanical instability of a helicopter is a complex phenomenon, involving cou- 
pling between the rotor and the body degrees of freedom. The rotor lead-lag regressive 
mode usually couples with the body pitch and/or body roll to cause the instability. 
When the helicopter is on the ground, such instability is referred to as ground res- 
onance. In the case of bearingless rotor, the ground resonance is usually avoided by 
providing elastomeric damper to the rotor blade and damping in the landing gear. Due 
to the inclusion of the elastomer, the dynamic characteristics of the bearingless rotor 
becomes very complicated. 
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Conventionally, the aeromechanical stability of a coupled rotor/fuselage system is 
analysed using a linearised perturbational technique. It is known that the elaistomer 
stiffness and damping characteristics are dependent on amplitude of motion. It is also 
highlighted in Sec. 2.1 that any linearised technique fails to capture these amplitude- 
dependent characteristics of the elastomer. The formulation and solution presented in 
this chapter aim to understjind the phenomenon of amplitude-dependent stability of the 
coupled rotor /fuselage dynamics under ground resonance condition. A methodology 
for analyzing the amplitude-dependent stability of the system under ground resonance 
condition is also presented. 


4.2 Formulation of Equations of Motion 

An idealised model of the coupled rotor/fuselage system is shown in Fig. 4.1. The blade 
is assumed to undergo only the lag motion. The elastomer is represented by non-linear 
spring and damping elements. The torque tube is assumed to be rigid and massless. 
The effect of pitch link in the lag mode is neglected. The hub is located at a height h 
above the C.G. of the fuselage. The fuselage is allowed to undergo rigid body pitch and 
roll motions. The characteristics of the fuselage support (landing gear) is represented 
by a linear spring and a linear damper. Since emphasis is placed on the effect of non- 
linearity of the elastomer on the ground resonance problem, the aerodynamic effects 
are not considered. The dynamical equations of motion of the coupled rotor/fuselage 
model consist of two sets of equations: one corresponding to the blade motion and the 
other representing the motion of the fuselage. 
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4.2.1 Blade Equations 

The equation of lag motion of an isolated blade is already presented in Sec. 3.2.2. Since 
the pitch and roll motions of the fuselage introduce additional inertia effects they must 
be included in the equation of lag motion of the blade. Accordingly, the kinetic energ>' 
of the k-ih blade in a A^-bladed rotor system can be written as 



where x (= 9yh) and y (= —6xh) are the perturbational displacements of the hub 
centre due to fuselage motion in pitch and roll modes, respectively (Fig. 4.1), 'ipk is the 
azimuth angle of the fc-th blade. 

The expression for strain energy due to bending is 



The strain energy of the non-linear spring due to the deformation A/ is 

= (4.3) 

where Keg is the equivalent stiffness which is dependent on the amplitude of motion of 
the blade. Hence the net strain energy of the system is 

U = Ui + U2. (4.4) 

The work done by the dissipative element of the elastomer is given by 

W'., = lc.,{A,)A," (4.5) 

where C^g is the equivalent (amplitude-dependent) damping coefficient. The equivalent 
stiffness and damping Ceg are taken as the dynamic stiffness {G') and damping 


80 



{G") for a given amplitude of the elastomer. Accordingly, from Eqs. 2.8 and 2.9, one 
can write 


G' = 
G" = 


ATe, = ATi - - ^K,X^ 


Ceq = 1 - 

U) 


4F 

ttujX 


where X is the amplitude of motion of the elastomer. 

The non-dimensional lag deformation of the blade is assumed to be of the form: 

^ T) = ^^i (^) ^ik {t) (4.6) 

i=l 

where is the zth rotating mode shape of the A:th blade, nl is the number of rotating 
modes. 

Since the equivalent stiffness of the elastomer is a function of the amplitude of motion, 
the rotating mode shape and frequencies will also be functions of amplitude of motion. 
The expression for 0, can be written as a sum of n assumed modes 


<i>i = ^iPlt + 02P2t + ••• + ^nPm 

where are the comparison function and pji is the t-th modal vector (J = 1, 2, 3...n) 
of the eigenvalue problem 

[if] {«} = A [M] {,} (4.7) 

with the elements of [K] and [K] given by 

and Mir = / ■ 

Jo 

The expression for ipl can be obtained following the procedure described in Sec. 3. 3. 2. 2. 
It is to be noted that since Keg depends on the amplitude, the eigenvalues and the 
eigenvectors are also dependent on the amplitude of motion of the elastomer. 
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Substituting Eqs. 4.1, 4.4 and 4.5 in Hamilton’s principle and neglecting higher 
order terms, the non-dimensional equation of motion of the ith lag mode of the ifcth 
blade can be written as 


nl nl ^ ^ r rl _ 1 

Zik {r) + Ztjt (t") + iijt - / UnOy sin V’i + ^hh cos = 0 (4.9) 

t=l r=l 


where C* = z = l,2...nl and k = l,2...N. 

Considering only the first lag mode {i = 1), the equation of motion of the m-th blade 
reduces to 


'Z\k ('^) + ^i^lk ('^) + 


6y sin 'tpk + cos 



(4.10) 


where A: = 1,2...A^. 

Applying the multi-blade-coordinate transformation, the above set of equations are 
converted into equations in rotor degrees of freedom. The collective and alternating 
modes are neglected as these do not participate in the ground resonance problem. The 
cyclic lag equations for the elastic lag mode can be wTitten as 


Zic ■+ Cz\c + 2iis + 

(uf - 1] 

1 ^Ic ^ ^Is P^x 

= 0 

(4.11) 

and zis -1- Czi 3 - 2iic + 


1 ^ ^ic P^y 

= 0 

(4.12) 

where P = ^h f ^id^ 

and 

C = 


(4.13) 


0 


4.2. 1.1 Fuselage Equations 


The in-plane (5„) and the radial (Sr) forces at the root of the A:th blade are given by 

Sy = J (-m^ + Tn^H'jdx- (i sin ipk-^ cos -ipk) mdx (4-14) 

and Sr = [ mQ^xdx - f 2Qtmdx -{£ cos rjjk + V sin “ipk) f mdx. (4.15) 
Jo Jo JO 
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The forces acting on the hub can be obtained by resolving the blade-root shear 
forces in the non-rotating frame and summing over all the //-blades. Considering the 
support structure as a combination of mass-spring-damper system having pitch and roll 
motions, the equations of motion of the support in, non-dimensional form, are written 
as 

+ Cy6y -t- KyOy ~ ^ 

and I A + CA + K^O:, - ^pzu = 0 (4.17) 

Zi 

Eqations 4.11, 4.12, 4.16 and 4.17 are arranged in the form 

[m] {1} + [c] {x} + [A-] {X} = 0 (4.18) 

where M, C and K are the mass, the damping and the stiffness matrices, respectively. 
The vector {X} represents the degrees of freedom, t.e., 

{X} = {z^c zu Oy . (4.19) 

The ground resonance stability is determined by the eigenvalues of the 8x8 matrix 
[G] where 

0 I 

[G] = . (4.20) 

4.3 Solution Procedure 

The step-by-step technique for solving the amplitude-dependent stability of the coupled 
rotor-fuselage system is described below: 

(i) For a given amplitude of motion of the elastomer, equivalent stiffness (ii'e?) and 
damping {Ceq) are evaluated. 
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(ii) Using Keg, the natural frequency {u>i) and the corrsponding mode shape (i^,) of the 
rotating blade in the lag mode are determined by solving the eigenvalue problem given 
by Eq. 4.7. 

(iii) Since only the first lag mode is considered in the ground resonance problem, using 
the mode shape the quantities C and p given by Eqs. 4.13 are evaluated. 

(iv) Knowing M, C and K (Eq. 4.18), an eigen analysis of matrix G is performed to 
determine the stability of the system under ground resonance condition. The eigenval- 
ues appear as complex conjugate pairs Sj — Oj ± iu)j where u)j represents the frequency 
of the mode and Oj represents the modal damping. The mode is stable if Cj is negative 
and unstable when Oj is positive. 

Since the lag mode shape is dependent on the operating conditions and the blade 
and elastomer parameters, for any change in these parameters the above steps have to 
be repeated. 

4.4 Results and Discussion 

A thorough investigation is carried out to study the influence of the elastomer on 
the ground resonance stability of a bearingless rotor. The elasomer is idealised as a 
combination of a non-linear spring, a Coulomb damper and a hysteretic damper. The 
parameters of the non-linear elastomer are given in Table 2.1. The data used in the 
analysis are given in Table 4.1. The values of the stiffness parameters of the support 
are chosen to provide frequencies in pitch and roll modes as 2 Hz and 4 Hz, respectively. 

4.4.1 Baseline Configuration 

For the baseline configuration, the locations of the elastomer and torque tube are set 
at = 0.1, and ^2 = 0.25. The damping in pitch and roll modes of the fuselage is 
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assumed to be zero (Cx = Cy — 0). The amplitude of motion of the elastomer is taken 
as 0.001 m. Following the steps mentioned in Sec. 4.3, the eigenvalues of the matrix 
G are calculated. Figures 4.2 and 4.3 show the variation of the modal frequencies 
(wj) and damping (aj), respectively, as a function of the rotational speed Q. Figure 

4.2 indicates that the coalescence of the regressive lag mode (zr) with fuselage modes 
in pitch (^j,) and roll (0x) occurs at three regions of rotational speed (namely 95-100 
r.p.m, 350-400 r.p.m and 440-600 r.p.m). In the range 350-400 r.p.m, 0y is unstable 
having maximum instability with cTj = 0.23 while in the r.p.m range 440-600, 6x 
is unstable with a maximum positive value of aj = 1.67 The coalescene of Zr 
with dy occuring at a very low^ r.p.m range 85-95 does not produce any instabilities in 
the system. The progressive lag mode {zp) is always stable with a slight variation in 
damping the levels. 

4.4.1. 1 Effect of Amplitude 

The ground resonance problem is solved for different three values of the amplitude of 
motion (0.0005 m, 0.001 m and 0.002 m) of the elastomer. The parameters are taken 
as Ci = Cy = 0, ^1 = 0.1 and ^2 = 0.25. The results, shown in Fig. 4.4, indicate that 
for this configuration, the damping in progressive lag mode decreases with the increase 
in amplitude. In other modes, the amplitude of motion of the elastomer has very little 
influence. 

4.4.2 Influence of Elastomer Location (^ 1 , ^ 2 ) 

A parametric study is undertaken to analyse the influence of and ^2 on the modal 
damping. The amplitude of motion of the elastomer is fixed at 0.001 m. 

The analysis is carried out by varying the value of keeping ^2 fixed at 0.25. Figure 
4.5 shows the results for four different values of ^1 (0.1, 0.07, 0.04 and 0.02). As ^1 
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maximum va.lue of damping coefficient C. These results indicate that by a judicial 
selection of the locations of the elastomer and torque tube, the stability of the system 
can be enhanced. 

4.4.2. 1 Effect Linear Support Damper 

The analysis is carried out for the baseline configuration (^i = 0.1, ^2 = 0.25). The 
damping ratios in the pitch and roll modes are assumed to be = .0529 and Qt — 
.082. The elastomer amplitude is taken as 0.001 m. Figure 4.10 shows that for the 
baseline configuration, the level of instability in the roll mode (0x) is a = 1.67 s“^ in 
the speed range of 440-600 r.p.m. With the inclusion of support damper, the maximum 
level of instability in roll mode is brought down to 0.87 s"^. But the region of instability 
spreads over a wider range of operating speeds (470-660 r.p.m). To highlight the 
importance of elastomer location on stability, the results corresponding to another 
configuration (^1 = 0.04, ^2 = 0.25) are also shown in Fig. 4.10. With support damping, 
this configuration has a maximum level of instability of 0.53 s~^ in roll mode which is 
lower than that of the baseline configuration. The region of instability is also reduced 
to 515-580 r.p.m. 

Keeping the other parameters same, when the support damping in the roll mode is 
further increased to = 0.20, the configuration = 0.04, ^2 = 0.25) becomes stable 
for the entire range of rotational speeds, whereas the configuration (^1 = 0.1,6 = 0-25) 
is not stable as shown in Fig. 4.11. 

Inclusion of the support damper in addition to the elastomeric lag damper reduces 
the level of the instability considerably. In addition, for a given value of fuselage 
damping, a proper choice of elastomer location can improve the stability of the system 
over the entire r.p.m range. 
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4.5 Summary 


An amplitude dependent stability analysis has benn carried out for a coupled rotr/fuselage 
system under ground resonance condition. It is observed that the elastomer amplitude 
does not alter the stability of the system in the rigion of maximum instabilities. How- 
ever, it exhibits appreciable influence on the progressive lag mode damping. The results 
also indicate that for a given location of the elastomer, there is an optimum location 
for the torque tube attachment from the point of view of ground resonance instability.. 
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Fig. 4.1 Idealised model of the coupled rotor/fuselage system. 
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Fig. 4.2 Variation of modal frequencies as a function of rotational speed. 
6 = 0.1, 6 = 0.25. 



Fig. 4.3 Variation of modal damping as a function of rotational speed. 
6 = 0.1, 6 = 0.25. 
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i (1/sec.) 



n (r.p.m) 

Fig. 4.4 Effect of amplitude of motion of the elastomer on modal damping. 

Amp. = 0.0005m; Amp. = 0.001m; Amp. = 0.002m. 

6=0.1,6 = 0.25. 
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Q (r.p.m) 

Fig. 4.5 Modal damping for different locations (^]) of the elastomer (zero support 
damping). 

6 = 0.1, ^2 = 0.25; 


6 = 0.07, ea = 0.25; 

6 = 0.04; 6 = 0.25; 

• ■ • • = 0.02, 6 = 0.25 
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'^'200 300 400 500 600 700 800 

Q (r.p.m) 

Fig. 4.6 Modal damping for different locations (^ 2 ) of the torque tube attachment 
(zero support damping). 

•••• 6=0.04.6 = 0.15; 

— 6 = 0.04, 6 = 0.23; 

6 = 0.04,6 = 0.35. 
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■ 200 300 400 500 600 700 800 

Q (r.p.m) 

Fig. 4.7 Modal damping for different locations (^ 2 ) of tho torque tube attachment 
(zero support damping). 

— 6 = 0 . 02 , 6 = 0 . 1 ; 

... ^1=0.02,6 = 0.2; 

.... 6 = 0.02, 6 = 0.3. 
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ty in roll mode for different combinations of 




(1/sec.) 



Fig. 4.10 Effect of support damping on stability = 0.082, = 0.0529). 

no support damping (^i = 0.1, ^2 = 0.25); 

with support damping (^1 = 0.1, ^2 = 0.25); 

with support damping (^1 = 0.04, ^2 = 0.25). 
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Table 4.1 Input data for ground resonance analysis. 

Parameters for non-dimensionalisation: 

mass per unit length, m = 9.7 Kg/m 

Length L = 6.6 m 

time Q = 32.8 rad/s 


Blade data: 

m 



9.7 Kg/m 

L 



6.6 m 

a 



0.0301 

Q 



32.8 rad/s 

Non- 

dimensional fuselage data: 


h 



1.0511 

h 



3.1023 

'Qr^ 

II 

h 

L 


0.2 



\l/2 

4 

Wx — 

U* 

/ 

Q 

U)y = 

(!<x 

y/2 

2 

n 

Ox 



0.082, 0.20 

Ov 



0.0529 

Cx = 

20x 
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Chapter 5 


TRANSIENT RESPONSE IN 
HOVER TO A STEP CONTROL 
PITCH INPUT 


5.1 Introduction 

Most of the aeroelastic response analyses of rotor blades address the problem of for- 
mulation and evaluation of the steady state response of the blade in forward flight 
for estimating the vibratory loads acting on the hub and the blade root. Very little 
information is available on the transient response characteristics of the rotor blades. In 
1953, Carpenter and Pridovitch [106] performed a wind tunnel experiment to examine 
the transient response beha^^our of a rotor blade, having only rigid flap motion. It was 
observed that for a step input in the pitch angle, the transient response of the blade 
in the flap mode exhibited dynamic overshoot leading to overstressing of the blade. In 
addition, the rotor thrust also showed a dynamic overshoot before reaching the steady 
state value. The dynamics of a bearingless rotor is much more complicated due to the 
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presence of non-linear elastomeric damper and multiple load path. The pitch input 
to a bearingless rotor blade, as opposed to that of a hingeless rotor blade, is provided 
by twisting the flexbeam through a vertical movement of the pitch control mechanism 
(Fig. 1.3a). Therefore, this configuration provides an additional kinematic constraint 
at the point P. It is expected that the non-linearity of the elastomer, multiple load 
path and kinematic constraint will have strong influence on the dynamic response char- 
acteristics of the blade. In addition, due to inherent differences in the construction of 
hingeless and bearingless rotors, there could be fundamental difference in their response 
characteristics. 

The aim of this study is to examine the transient response to a step input of an 
isolated bearingless rotor blade in hovering condition simulating a wind tunnel test. 
An attempt is also made to quantify a measure of aeroelastic couplings (pitch-flap 
and pitch-lag couplings) from the steady state response and the dependence of these 
couplings on various geometric parameters of the bearingless rotor blade. 

5.2 Formulation of Equation of Motion 

The transient response analysis of a bearingless rotor blade requires the formulation 
of aeroelastic equations of motion representing coupled axial, flap, lag and torsional 
motions of the blade, including non-linear elastomer and kinematic constraint. The 
equations of motion are derived using Hamilton’s principle. The important ingredients 
required for the formulation are: 

(i) the strain energy of the blade, 

(ii) the kinetic energy of the blade, 

(iii) the strain energy and dissipative energy associated with the elastomer, 

(iv) the work done due to aerodynamic loads 


100 



and 

(v) a suitable model of the kinematic constraint. 

5.2.1 Assumptions 

Since the aim of the study is to understand the physics of the dynamical motion 
and aeroelastic couplings in a bearingless rotor due to multiple load paths, pitch link 
geometric constraint and non-linear elastomer, several assumptions have been made 
while formulating the equations of motion. These are: 

(i) the rotor blade is an Euler-Bernoulli beam made of homogeneous isotropic ma- 
terial, and undergoing coupled axial, flap, lag and torsional deformations. 

(ii) the torque tube and pitch link are rigid and massless elements (usually, the 
torque tube has a very high stiffness relative to that of the flexbeam or rotor blade 
section), 

(iii) the rotor blade has zero precone, 

(iv) the cross-section of the blade is symmetric, 

(v) the pitch control mechanism is infinitely stiff, 

(vi) the rotor is operating under hovering condition simulating a wind tunnel test, 
and (vii) the aerodynamic loads on the rotor blade are represented by a quasi-steady 

two dimensional aerodynamic model with a time varying inflow. 

The mathematical procedure of formulating the equations of motion essentially 
follows the general steps adopted in Refs. [35, 46]. 

5.2.2 Ordering Scheme 

In the formulation of equations of motion of a rotor blade, a large number of higher 
order terms are generated. In order to identify and eliminate the higher order terms in 
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a consistent manner, an ordering scheme is employed. This ordering scheme is based on 
the assumption that the slope of the deformed elastic blade is moderate, and of order 
e (0.10 < e < 0.20). Orders of magnitude in terms of e are then assigned to various 
non-dimensional parameters governing the rotor blade problem. In this formulation, 
terms of order ^ are neglected with respect to terms of order 1, t.e., 

0(1) + 0(e2) «0(1) 


The order of magnitude of various non-dimensional parameters governing this problem 
are given below': 


ORDER 1 

I r_^ ^ _ 1 a. 
i’ ^dx' dr ~ ildt 

ORDER 

ORDER e 


n. L v yi y -in 6 6 0 9^ 
ORDER 6^ 


u ^ 

I, u.,X) txn 


All symbols are explained in the nomenclature included at the beginning of the thesis. 


5.2.3 Coordinate Systems 

The description of complex deformation of a rotor blade requires several coordinate 
systems. Each coordinate system is represented by a triad of orthogonal unit vectors. 
Coordinate systems used in deriving the equations of motion are described below: 


102 



5. 2. 3.1 Hub Fixed Rotating Coordinate System 

The {ej,ey,ez) system shown in Fig. 5.1 rotates about Z-axis with constant angular 
speed Q of the rotor. Its origin is fixed at the hub centre O and the direction of 4 is 
along the elastic axis of the blade. 

5. 2.3.2 Undeformed Cross-section Fixed Rotating Coordinate System 

The (eiijCyi, 621) system is parallel to the hub fixed rotating coordinate system. But 
its origin is located at O' at a distance x from the hub centre O as shown in Fig. 5.2. 

5. 2.3.3 Undeformed Curvilinear Cross-sectional Coordinate System 

The curvilinear {e^, ijj, e^) system has unit vectors e,, and along the principle axes of 
the blade cross-section. The orientation of (e,,,e^) with respect to (6^1,621) or (ej„e2) 
is represented by an angle 0 (Fig. 5.3), where 6 is the sum of pretwist and pitch 
input 6,n. 

5. 2.3.4 Deformed Cross-sectional Coordinate System 

The (ex', Cy', 62') system represents the deformed cross-sectional coordinate system. It 
is obtained after performing a rigid body translation and rotation of (eii,'eyi, 621), 
corresponding to the local axial, flap and lag deformations at the elastic axis. The 
direction of ix' is along the tangent to the deformed elastic axis. This system is used 
to define the directions of flow velocities for aerodynamic load calculations (Fig. 5.4). 

5.2. 3.5 Deformed Curvilinear Cross-sectional Coordinate System 

The {e^' , Cjj' , system represents the orientation of cross-sectional principal coordi- 
nates after axial, flap, lag and torsional deformations. The angle between the coordi- 


103 



nate systems (e^', e^/, ey) and represents the sum of pretwist (0tw), pitch 

input (^,„) and elastic twist of the blade ((^), as shown in Fig. 5 . 5 . 

5. 2. 3. 6 Local cartesian Coordinate System 

The system (j/i,y2)?/3) represents a local cartesian coordinate system with its unit 
vectors (61,62,63) parallel to the orthogonal triad (ei,e,„ej of the cross-section. The 
Green strain measures are obtained in the undeformed curvilinear coordinate system 
(ix, e,j, 6^). These have to be transformed to the local cartesian system because the 
constitutive relation of the material is defined in the local cartesian system. The 
essential difference between {ex,e^,e^) and (61,62,63) is that the derivatives of e,,, 
with respect to x represents the twist rate of the blade in the undeformed state, whereas, 
the derivatives of 62, 63 with respect to x are both 2;ero. 

5.2.4 Blade Model 

The strain energy and the kinetic energy expressions of the blade are detailed in the 
following sections. 

5.2.4. 1 Strain Energy of the Blade 

Details of the derivation of strain components are described in Appendix A. Consistent 
with the assumption of Euler-Bernoulli beam, the non-zero strain components at any 
arbitrary point (77, C,) on the cross-section of the blade can be written in terms of u, v, 
w and <j> as 

Cxx = «x + Y + + + 

- V^xx [v COS {9tw + 9,n + <!>)-( sin {$tw + Oxn + <l>)] 

- w^xx [v sin {9tv, + + <?i) + C cos {9tu, + + 4^)] ' (o-l) 
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(C “b ^»?) 0,1 

(5.2) 

■'xC = (77-^c)0,x 

(5.3) 


where u, v and u' are the elastic displacements of any arbitrary point on the elastic 
axis in the Cx, e.y and directions, respectively. The elastic twist is (j> and ^ is the 
warping function for the cross-section due to torsion. 

The strain eneg^' of the blade is given as 


^--\l IL i^xx^xx + OxTilxT, + cTxO^d dVOL 


(5.4) 


where 


dVOL = ^/gdr]dQdx 

XX — ^^xx 
O'xjj = G'^xrj 

CTxC = G-IxQ 

and g is the determinant (equal to 1) of the metric tensor defined in Eq. A. 28 in 
Appendix A. 

The variation of the strain energy is given by 

6U = i:iL i^^xx^^xx + CTxnS'TxT, + (^xc^hxd dr}dC,dx (5.5) 

where 

5txx - du^x + v^xdv^x + W,xdw^x + + C^) {^tw + ^in + <t>) ~ '^^(p^xx 

- [t] cos {9tyj + 9in + (p)-C sin (^^ + + <P)] {dv,xx + ‘W,xxd<t>) 

\j} cos (^Oiyj Bin -|" 0) -f" ^ sin {Olyj -H Om "f" 0)] i^dw^xx ^,ix^0) 

hxv = -C<^0,x 
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fj = 77 -^^ 

C = C + 

and the subscript , denotes a partial derivative with respect to the subsequent sub- 
script (s). 

Performing cross-sectional integration and applying the ordering scheme, the vari- 
ation of strain energy expression can be written as 

rL 

SU / {y^x “t" “I" ^0 P e'^^,xx 

J 0 

+ \Mz> cos {dt^ + 6i^ + <^) + Myf sin {6tw + + o)] {Sv^xx + (5.6) 

*4" sin *4” ^in 4“ cos {Ofxjj *4“ 4“ 0 )] (^^,xx ^,xx^^)} 

where 

Vx' = EA -I- -f 

P [^,xx COS {0-iyj 4“ ^tn 4" 4* w^xx sin 4" Olfi "4“ ^)] 

Sx' = GJ(j>^x 

( v\ w\ \ 

Tx^ = EACo4>.xUx + ^ + -f\ 

Px' “ P ^Pz^yXx 4* PAD2 [*^,xx COS (Ofiij 4“ ^tn) ^ ,xx i^tw 4” ^in)] 

Myt = Plfpij [t^^XX (^tiy 4“ 6ifi 4" ^,XX (^ft£' 4~ Oiji -f- ^)] J5yl£?2^^,XX 

and Mz' = P^C [^»xx cos {9tw 4- Oin 4- </>) 4- w^xx sin {Ot^- + 4- (f))] 

/ y2 yj2 \ 

- EA%L,+ f + -f\-EAQ{e,.. + e,„)^<,^ (6.7) 

The sectional integrals in Eqs.(5.7) are defined as follows: 

EA = J J EdrjdQ 

EArja = f J [Ev]dvdC 
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El„ = j f [£C"] dT]d( 

EIk = j j 

EACo = j J [E (rf + C^)] dT]d(^ 

GJ = 1 1 [G {f}^ + C^)] dvdC 
EAC 3 = J J (t]^ + C^y^dT]d(; 

EACi J J [^V + C^)] drjdC 
EADz = f J dvdC 
and EAD 2 = II [£;C^] drjdC 

Because of the antisymmetric character of ^ and the assumed symmetry of the cross- 
section about 77 -axis 

f J [£4']<i7)<iC = 0 

f f [EQdndi^^ 

I j [BC (>;=' + C^)] dr,d^ = 0 
j j [Bi' {rf + c")] dr,dl; = 0 

and j j[E7]^]drid(; = 0 


5. 2.4.2 Kinetic Energy of the Blade 

Combining Eqs. A. 17, A. 14 and A. 16 in Appendix A, the position vector of any point 
P on the deformed beam is given by 


R = (x + u) Ci + viy + wcz + Tjejji -f Ce^' + a^ex> 


= + Vl^y + ZlBz 


(5.8) 
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where a is the warping amplitude. Noting that the warping amplitude a is taken as 
equal to elastic twist rate of the blade (Eq. A.35) and using the transformation 
matrices, the expressions for xi, yi and zi can be written as 

Xi = a: + u - u_,r[7?cos(0tu, + + (^) - Csin(0tu, + + (;^)] 

- w^x [rj sin {9trv + 0,n + ^) + C cos + 0i„ + (p)] 
yi = v + r]cos{etnj + Oirt + (f))-Csm{ety, + ein + (i)) (5.9) 

and zi = m + r? sin {Ot^j + 9in + <f>) + ( cos {dtru + Oin + 4>) 

or, 

xi = X + u - - v^x ivi -v)- W,x (Zi - w) 

yi = t; + (yi-u) (5.10) 

and zi = w + (zi— w) 

where the terms (j/i — v) and {zi — w) are introduced for convenience. 

The velocity of the point P is given by 

V = ^ + nexXR 
at 

= (xi -flt/i)ei + (yi4-flxi)ey + iiez (5.11) 

where 

xi = u - [yi -v)- (w,x - i) (zi -w)- '^<p,x 

yi = v-^{zi-w) 

and ii = w + ^{y\-v) (5.12) 

The kinetic energy of the beam is given by 

r(f) = ^ /V I 
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The variation of kinetic energy is 


5T = I IL pV6Vdr)d(^dx ( 5 - 14 ) 

Substituting for V from Ekj. 5.11 and integrating 5T by parts with respect to time, 
one gets 

ST = f [ [ P + ZySv + Z^5w + Z^’Sv^i + Z^>6w^x + Z^6(f)] dT]d(^dx (5.15) 

0 d J A. 

where the terms Zu, Z^, Z^, Z^t, Z^' and Z^ Eire the coefficients of Su, Sv, 6w, Sv^x, 
5w^x and respesctively. 

Integration of the expression over the cross-section yields 

= ^ ^Z^Su + ZySv + Z^^Sw + Z^'6v^ + Z^iSw^x + Ztj)5(j)^ dx ( 5 . 16 ) 

After eliminating the higher order terms and assuming the cross-sectional symmetry 
about 77 -axis and an antisymmetric warping function the expressions for Zu, Z^, Zu,, 

Z^> , Z^' and Z^ can be written as 

Zu = m (Q^x -I- 2 f 2 uj 

Zu = -f rriTjm cos {Otu, + 0^ + 4>)- 2mQu 

+ 2rm]m^ [v^x cos {9tuj + 0 ,„) -f sin (^^ + ^in)] - mv + rnrjm^ sin (dtw + ^m) 
Zyj = -mw - 7n77nj^COS {9tui + 9in) 

Z^ = Irnrf,^)cos{dtu~-\r9,n + <i>)sm{dtui + 9m + 4)) ( 5 . 17 ) 

- mrim^'^x [w,x cos {9tu, + ^,„) - sin {9tw + 9^)] - mrjm^'^v sin {9tu- + 9,n) 

- mr)m [V sin {9tw + ^in) - W cos (0tu. + ^in)] 

Zu' = cos {9tw + 0,n + '^) + cos {9tui + + <^)] 

and Z ” ^ 277^71 X sin i^9^u) “h ^) d" 2^v sin (^9 tw + 9in + ^)] 
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where sectional integrals are defined as 


m 

^ 1 1 

mrim 

= f J Pd'ndQ 


= j J f<^dT]d(: 


= J J PV^dTjdC 

Im 

II 

+ 

>->< 

3 


The integrals involving the terms (yi — v) and {zi - w), introduced for convenience in 
Eqs. 5.10, are given by 

J j p{yi-‘^) dvdC = mrim cos + <f>) 

J j p{zi-w) drjdC = mrjm sin + 6in + (f) 

J j p{yi-v) (zi - w) drjdC = - Im^ cos (Ot^ + Sin + 4>) sin + 6in + (f>) 


J j P [(yi - vf + (-21 - ic)^] = dm 


5.2.5 Elastomeric Damper Model 

An elastomer is placed between the torque tube and the flexbeam to provide adequate 
lag damping to the rotor blade. It also serves the purpose of a spacer between the torque 
tube and the flexbeam. During blade deformation, the elastomer provides constraints 
in flap, lag and torsional motions. The elastomer data are usually available [25, 26] 
for the lag mode only. Therefore, following the approach of Ormiston et al. [31]^ the 
elastomer is represented by a very stiff linear spring in the flap mode. In the lag mode, 
the elastomer is modelled as a parallel combination of a non-linear spring, a Coulomb 
damper and hysteretic damper (Fig. 2.6); and the effect of elastomer is not considered 
in the torsional mode. 
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5.2.5.1 Strain Energy of the Elastomer 

Flap Mode : 

The expression for the strain energy of the spring due to deformation A/ is 

U, = iji-.A) (5.18) 

where A/ = w (L 2 , t)-(L 2 -Li) - m (Li,t) 

Variation in spring energy is given by 

5Uf = KiUfShf (5.19) 

Lag Mode : 

The expression for the strain energy of the spring due to deformation A/ is 

U, = iifiAf-ijfaAf+lKsAf-iA'rAf (5.20) 

where Ai = u(L 2 ,<) - (ij - ii) - u (£'].*) 

Variation in spring energy is given by 

5Ui = (iiTi A, - + JTsAf - JTtAJ) <5A, (5.21) 

5. 2. 5.2 Work Done by the Dissipative Element 

The expression for the virtual work done by the dissipative (Coulomb and hysteretic) 
elements of the elastomeric damper in the lag mode is given by 

8Wi = Frf<5A( (5.22) 

where Fd = —A/ + Fsgn | A; | 
u 
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5.2.6 Aerodynamic Model 

The aerodynamic forces and moments are treated as external loads acting on the blade. 
The expressions for aerodynamic loads are not obtained explicitly, rather they are 
derived by an implicit formulation, following the procedure adopted in Ref. [52]. The 
aerodynamic loads are included separately in the computer programme. At each time 
step, these loads are evaluated knowing the response of the blade and other velocity 
parameters. 


5.2.6. 1 Aerodynamic Loads 


The time varying aerodynamic lift and pitching moment acting on a typical cross- 
section of the blade are evaluated based on Greenberg’s extension of Theodorsen’s 
theory [112] developed for a two-dimensional airfoil undergoing harmonic motion in a 
pulsating incompressible flow. A quasi-steady approximation of the unsteady theory is 
used, by assuming the lift deficiency function c{k) = 1. 

Based on the quasi-steady approximation [112], the expressions for aerodynamic 
lift La and pitching moment Ma per unit span are given by 


-aipAbl 


d {hf + YfO) 


dt 




-t- aiPAh,Vf[h}^Vfe+{hx-XA)e] 


(5.23) 




bi 


+ 


a.pAbiXAVf [hf -f Vfd + (6i - Xa)( 


(5.24) 


where oi is the two dimensional airfoil lift curv^e slope; bi is the semi-chord ; pa is the 
density of air; hf is the plunging velocity; Vf is the free-stream velocity; Xa is the 
blade offset of the aerodynamic centre from the elastic axis. The pitch angle 0 with 
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respect to the free-stream is given bj' 



Substituting Eqs. (5.28) and (5.29) into Eqs. (5.23) and (5.24), the expressions for 
aerodynamic lift and pitching moment per unit span are written as 

La = -Uc-(xA-^y 

+ aiPAb,{V„.-Uf.e)[-U^+{bi-XA)9] (5.30) 

Ma = la,pAblUxA-^j^u^' + lb,{u,.-u^.e)e+ «j 

+ aiPAbiXA{U^--U^O)[-U^,+(b^-XA)il] (5.31) 

The expression for the profile drag per unit span is given by 

D = Ca^PaKVI = CA,PAh (t/J + V}) (5.32) 

Evaluation of the sectional hft and moment require U^’, U^, 6 and 6 which are 
obtained first as outlined below. 
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Blade Velocity Relative to Air 


The velocity vector of a point on the elastic axis of the blade relative to air is given 

by 

U = Vo-Va 

— Ux’Cx' + Ufferji + U^ie^i ( 5 . 33 ) 

where Vo is the velocity vector of a point on the elastic axis of the blade and Va is 
the velocity vector of air due to inflow only, since the present analysis is performed for 
hover. 

The velocity vector of a point on the elastic axis of the blade, Vo, can be obtained 
from Eq. 5.11 as 


Vo = V,=c=o = + V^ty + y/e, 

where 


[ K’ ' 


G 

1 

•3 

\ 

I/O 

y 

► = i 

fix + {? + fiu ^ 

K" 

V ^ ^ 


w 


The expression for Va can be given by 

VA = V*e, + V*e, + V,‘‘e. 

where 


V/ 


f \ 

0 

V/ 

> = <( 

0 ► 

\ ^ J 


—V 

\ / 


( 5 . 34 ) 


( 5 . 35 ) 


( 5 . 36 ) 


( 5 . 37 ) 
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Where v is the inflow velocity through the rotor disk. The velocity components U^' 
and U(-i are obtained by combining Eqs. (5.33), (5.34) and (5.36) as 


Ux' 


' v° - v-/ ■ 

Ur,' 

. = [T] . 

yo _ t/a 

V y 



y/ - V/ 

^ z * / 


(5.38) 


where the transformation matrix, [T], between the deformed curvilinear co-ordinate 
system {e^' ,e ^' and the hub fixed rotating co-ordinate system {ex,ey,ez) under 
small angle approximation, is given by Eq. A. 12. 

The components of the aerodynamic forces and moments per unit span in the 
deformed co-ordinate system are related to L and M (Fig. 5.5) by 


Pr,' 


= LsinaA — D cos a a 
p^t = LcosaA + D sin a a 
9x' = ^ 

where a a is the local angle of attack of the blade and given by 


=r —tan 


-1 


'U, 


Sin a A — 


cos a A 


cl - 


u. 


c' 


Ur) 


'U, 


!L - 




(5.39) 

(5.40) 

(5.41) 


(5.42) 

(5.43) 

(5.44) 


\Ur) y/u/Tup 

The aerodynamic forces and moments per unit span in the hub fixed rotating co- 
ordinate system (e*, 6^,6^) can be obtained as 


(5.45) 


' N 

Lu 

1 Lv 

. = [Tf < 

1”! 

Pr}' 



, Pc' - 
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and 


Qx 


Qx' 



0 



0 


(5.46) 


The effect of the non-conserv'ative distributed forces and moments due to the aero- 
dynamic loads are included by applying the principle of virtual work. In general, the 
components Qy and Qz are of much smaller order compared to and Hence 

they are neglected in the formulation of virtual work. 

The virtual work done by the distributed aerodynamic forces is given by 


{SWe)f = J (LySv 4- L^Sw) dx 


(5.47) 


Similarly, the virtual work done due to the distributed moment is given by 

f%Mdx (5.48) 

J 0 

Thus, the total virtual work {5We) due to aerodynamic forces and moments can be 
written as 


sw, = + {SW,)m 


(5.49) 


5.2.6.2 Inflow Model 

Since the present study addresses the transient response of a rotor blade in hover, 
it is essential to model the time varying nature of inflow through rotor disk. In the 
quasi-steady aerodynamic model presented above, the rotor inflow (i/) is included as 
a parameter which has to be supplied at every instant during the transient analysis. 
It is known that the blade motion causes unsteady wake, w'hich in turn affects the 
blade loading. The dynamic inflow model [113, 114] simulates the unsteady wake 
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effects in a simple form. In this model, the inflow variation in collective mode (all 
blades have identical response) is assumed to be uniform over the entire rotor disk. 
The dynamic inflow model has been successfully applied to several aeroelastic and 
aeromechanical stability problems [83, 85, 115]. An extension of the dynamic inflow 
model is the generalised wake model [116]. The generalised wake model takes into 
account the azimuthal and radial variations in the inflow velocities. Manjunath et 
al. [117, 118] applied the generalised wake model to analyse the flap-lag and flap-lag- 
torsional stability problems in hover and forward flight. 

In the present study, a simple dynamic inflow model has been used to evaluate 
the inflow at every instant of time from the knowledge of the blade loading. In the 
present problem, a step input is given to an isolated blade in hover simulating a wind 
tunnel test, which essentially corresponds to the collective mode operation of the rotor. 
(It may be noted that in collective mode, all the blades in the rotor system execute 
identical motion). In the collective mode condition, the relation between the rotor 
thrust and inflow is given by [85, 113, 119] 

T = Thu + 2'kL^pai^ (5.50) 

where m = 0.637/9A^trT^ (apparent mass) 

and u = induced inflow velocity uniform over the rotor disk. 

Assuming there are N blades in the rotor system, the expression for T is also given by 

f = f Y^(U),dx (5.51) 

.=1 

where the subscript i represents the ith blade. Since in the collective mode in hover, 
all blades experience identical loading, one can write 

T = / NLyjdx , (5.52) 

Jo 
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where is given by Eq. 5.45. 


The non-dimensional thrust coefficient Ct is given by 


Ct — 


f 

PATrL^iClL)^ 


(5.53) 


5.2.7 Kinematic ConstrEiint due to Pitch Link Motion 


Figure 5.6 shows a bearingless rotor blade with pitch link at the leading edge. For a step 
pitch control input, the point P on the pitch link is moved by the pilot in the vertical 
(e^) direction through a distance Wp and is held at that position. Since it is assumed 
that the pitch control system is rigid, the point P cannot have any displacement in 
the vertical direction (e^) during dynamical motion of the blade. It is to be noted that 
the point P is not constrained in other five degrees of freedom, namely displacements 
along Ci, Cy and rotations about Ci, Cy, e^. From kinematics, the vertical displacement 
at the point P due to the blade motion can be given by 


Wp 

where wp 
Wc 
w'c 
and 0c 


Wc — w'cb -I- 0 co 

vertical displacement of the pitch link, 
flap bending deflection at the point C, 
slope at the point C 
torsional deflection at the point C. 


(5.54) 


The imposition of the kinematic constraint corresponds to satisfying the condition 


Wp = Wp 


(5.55) 


at all instants of time. 

Due to this kinematic constraint, a reaction force Fp in the vertical direction is 
developed at the point P. In the response analysis, this force is transferred to the 
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blade as an external loading at the point C. The external loads (bending moment , 
a torsional moment Mf and a concentrated force Ff) are given by 



H 

II 

(5.56) 


= —FpbCy 

(5.57) 

Mf 

= FpaCx 

(5.58) 


where ‘a’ is the pitch horn chordvdse offset and ‘6’ is the distance of the pitch link from 
the section AA along the length of the blade (Fig. 5.6). 

The values of Fp, and are changing continuously so as to maintain the 
constraint mentioned in Eq. 5.55. The virtual work (SWunk) by Fp, Mf and are 
given as 

SWunk = F^Su'c + M^Swc^x + Mf6<j)c (5.59) 

5.2.8 Treatment of the Axial Degree of Freedom 

In solving aeroelastic stability and response problems of rotor blades, the treatment of 
axial degree of freedom of the blade is ein important consideration in order to properly 
account for the centrifugal and coriolis effects associated with the axial mode. In the 
literature, there are three different approaches, namely, (i) elimination of the axial 
degree of freedom, (ii) substitution technique and (iii) use of axial inextensionality 
criterion. The details of these approaches are given in Refs. [58, 120]. 

In the present formulation, the axial inextensionality criterion is used to eliminate 
the axial degrees of freedom from the equations of motion. In this technique, the axial 
strain at the elastic axis is set equal to zero. This is based on the physical consideration 
that the blade is axially very stiff. Thus, the axial strain at the elastic axis (Eq. A. 32) 

= + + = 0 ■ (5.60) 
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With this assumption, the axial displacement u at any distance x from the root of the 
blade, is expressed in terms of the bending displacements (v and w) as: 

“ = + (5-61) 

The time derivative of u is expressed as 

u z= - f 4- w,:cW,x) dx (5.62) 

J 0 

Using Eqs. 5.61 and 5.62, u and its time derivative are expressed in terms of the 
bending deformations in the energy expressions. By this process, the axial degree of 
freedom is completely eliminated from the equations of motion. 


5.2.9 Equations of Motion 


Using the variations of energy expressions given by Eqs. 5.6, 5.16, 5.19, 5.21, 5.22, 5.49 
and 5.59, Hamilton’s principle can be written as 



-6U- SUi - SUf + 5Wa + + SWunk) dt = 0 


(5.63) 


Following the assumed mode method, the flap, lag and torsional deformations of 
the blade are expressed in terms of the generalised coordinates (r,, 9,) and shape 

functions (<^i, <^i, ji) as 

i=l 

n2 

i=l 

n3 

i=l 

X 

where r = fit ,(=y 
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The equations of motion are derived by substituting Ekjs. 5.64 in the energy expression 
and taking variations of the generalised coordinates. These equations are non-linear, 
coupled ordinary differential equations which can be expressed in symbolic form as 

[M] {y} + [C] {y} + IK] {y} + (y,y)} = {q (r,y,y)} (5.65) 

where M , C and K are the mass, damping and stiffness matrices, respectively. All the 
non-linear terms are included in the vector Q, represents the generalised force 
associated with the aerodyncunic loading. 

5.3 Solution Procedure 

The transient response of an isolated bearingless rotor blade to a step control input 
in pitch angle is obtained by simultaneously solving the non-linear coupled equations 
(Eqs. 5.65) and the inflow equation (Eq. 5.50) by numerical integration satisfying the 
kinematic constraint given in Eq. 5.55 at every time step. 

The step-by-step procedure of the numerical technique is given below: 

(i) At t = 0, the point P of the pitch link (Fig. 5.6) is given a vertical displacement 

Wp. 

(ii) An initial value of the pitch link load Fp (t) is assumed. Knowing Fp, the external 
loads {Fp, and Mf) acting at the junction of the flexbeam and blade are calculated 
using Eqs. 5.56, 5.57 and 5.58. 

(iii) At each time step, knowing the value of Fp{t), blade response and inflow, the 
following sub-iteration is performed to estimate the correct value of Fp{t+At) satisfying 
the kinematic constraint, at time step t + At. 

(a) The inflow velocity and blade response are calculated at t -I- At. 

(b) The elastic deformations at the point C {wc, wc,x 2 uid (pc) are evaluated and 
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using Eq. 5.54, the vertical displacement of the pitch link Wp (corresponding to Fp) is 
calculated. 

(c) Fp is incremented by an amount AFp. For this value of Fp + AFp at time t, inflow 
velocity and blade response are obtained for time t + At. 

(d) Knowing the blade response, the modified displacement, wf (for Fp + AFp), of 
the pitch link is evaluated. 

(e) The proper value of Fp{t + At) is estimated by linear interpolation of w$ (for Fp) 
and TUp (for Fp + AFp) to match Wp. 

(iv) With the new value of Fp{t + At) and the response of the blade and inflow at 
t + At, step (iii) is repeated until a steady state is reached. 

It may be noted that the response of the blade and inflow at t + At correspond to the 
constraint load Fp{t) at time t. 

To initiate the response calculation starting from time t = 0, an initial estimate 
of Fp{0) is required. This value of F|>(0) is estimated from static analysis considering 
isolated flap bending and twisting of only the flexbeam portion of the blade. The out 
of plane deflection of the flexbeam at the point C due to a concentrated load Fp and 
bending moment applied at C is taken to be that of a tip-loaded cantilever, given 
by 

FpLl M9LI 
- ,.. T ^ q e — £ 

^ ZEI 2EI 


Substituting Mf from Eq. 5.57, one obtains 

FpLl Fpbhl 


Wc 


(5.66) 


3EI 2EI 

where L 2 and El are the length and flexural rigidity of the flexbeam. 

The slope of the flexbeam at the point C due to a concentrated load Fp and bending 
moment is given by 

FpLj FpbLt 


Wc 




2EI 


El 


(5.67) 
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(5.68) 


The elastic twist (j)^ at the point C due to the twisting monient is given by 

I Ij2 ^ FpciL2 

“ GJ " ~gT 

where GJ is the torsional rigidity of the flexbeam. 

Combining Eqs. 5.66, 5.67 and 5.68, the following relation between Wp and Fp is 
obtained 


FpbL\\ ( FpL2 FpbL2\ , FpaL2 
^ V ZEI 2EI ) 2EI ~b7 " ) 


(5.69) 


From Eq. 5.69, an initial estimate of Fp can be obtained for a given value of w‘p. 


5.4 Results and Discussion 

The transient response analyses are carried out for two types of rotor blade configu- 
rations, namely, (i) hingeless blade and (ii) bearingless blade. In both cases, a step 
pitch input is applied to the blade operating in the hover condition simulating a wind 
tunnel test. The transient response analysis of a hingeless blade is carried out with the 
following objectives: 

(i) to validate the equations of motion and the solution procedure 
and 

(ii) to compare qualitatively the dynamic response characteristics of hingeless and 
bearingless rotor blades. Such a comparison will provide a usefull insight about the 
behaviour of these two different rotor configurations. 

5.4.1 Hingeless Rotor Blade 

For the hingeless blade configuration, the formulation of equations of motion essentially 
follows the same procedure outlined in Sec. 5.2. The terms involving th6 strain energy 
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of the elastomer (Eqs. 5.18 a.nd 5.20) and the work done by the damping element of the 
elastomer (Eq. 5.22) are not included in the Hamilton’s principle (Eq. 5.63). Since the 
pitching motion to the blade is given at the root (Fig. 1.2a), no kinematic constraint 
(Eq. 5.54) is imposed on the blade. 

The assumed mode shapes (Eqs. 5.64) for bending and torsional deflections are the 
undamped and uncoupled free vibration modes of the rotating blade. In the response 
analysis, the six modes, namely, three modes in flap, two modes in lag and one mode 
in torsion are considered. The number of blades in the rotor system is taken as N = 4. 
The data used for numerical calculations are given in Table 5.1. 

A step input of 0.1 rad (0i„ = 0.1 rad) is provided at the blade root. The blade and 
inflow equations are solved simultaneously by a numerical integration routine (NAG) 
based on Runga-Kutta-Marson algorithm. To check the validity of the time step, time 
integration was carried out with three time steps viz.. At = 0.0015 s. At = 0.0035 s 
and At — 0.006 s (in real time). In all the three cases, the results are found to converge 
to the same values. In all the subsequent computations, a time step of At = 0.003 
s is used. The response calculations are continued till the steady state condition is 
reached. The results of the transient response are shown in Figs. 5.7a-c. In Fig. 
5.7a, the non-dimensional flap response at the blade tip (wtip) is plotted against non- 
dimensional time (r). An overshoot in the flap response is observed the initial phase 
of the response. Subsequently, the flap response reaches the steady state value 0.019. 
The maximum value of overshoot is found to be (0.055), which is approximately 2.3 
times the steady state value. This kind of overshoot in the transient response in the 
flap mode of a rotor blade was also observed in the experimental results reported in 
Ref. [106]. 

The variation of the lag tip response {vup) with r is shown in Fig. 5.7b. Since 
aerodynamic damping in the lag mode is small, the response exhibits a large amplitude 
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of oscillation with a very low decay rate. Figure 5.7c shows the variation of torsional 
response {(pup) at the blade tip. The steady state value of the elastic twist is found to 
be -0.005 rad. with a very small overshoot during the initial phase. The variation of 
inflow velocity shown in Fig. 5.8 indicates that it starts from zero and reaches a steadv 
state value 0.048 with a little overshoot. The thrust generated by the rotor, shown in 
Fig. 5.9, behaves like the flap response. The maximum value of thrust coefficient is 
0.0057 which is 2.28 times the steady state value. 

The variation in the maximum thrust overshoot with respect to the rate of pitch 
input has also been studied. The transient response calculations are performed with 
four different rates of blade pitch input. The initial rate of pitch inputs are 20 deg/s, 
48 deg/s, 120 deg/s and 200 deg/s; after reaching a value of 0.1 rad. the pitch input 
is held constant. In each case, the ratio of the maximum to the steady state 

values of rotor thrust (Ct,) is obtained. Figure 5.10 shows the variation of 
Vs. rate of pitch input. It is seen that initially the ratio increases with 

the blade pitch rate and reaches a steady state value 2.37 at high pitch rates. 

Carpenter and Pridovich [106] carried out experiments on the transient response 
of a rigid rotor blade with varying pitch rate inputs. The blade is allowed to undergo 
only flap motion. They correlated their analytical results to the experimental results. 
(It may be noted that in Ref. [106], Carpenter et al. have not provided any blade 
data). Their experimental results are shown in Fig. 5.11 where thrust overshoot for a 
blade pitch rate of 200 deg/s is about 2 times the steady state value. In the present 
analysis, for the same blade pitch rate (200 deg/s), the overshoot is found to be about 
2.37 times the steady state value. This difference in magnitude may be attributed 
to the differences in the blade model. The qualitative nature of the present results 
conforms with the experimental observation, thereby validating the blade equations, 
inflow' model and the solution procedure. 
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5.4.2 Bearingless Rotor Blade 

A thorough investigation is carried out to study the influence of the non-linear elas- 
tomer and kinematic constraint on the transient response characteristics of an isolated 
bearingless rotor blade to a step input in hover condition simulating a wind tunnel test. 
The parameters of the elastomer are given in Table 2.1. Since the transient response 
of the bearingless rotor blade is compared with that of the hingeless blade, the same 
blade data is used. It is important to point out that in a practical bearingless rotor 
blade configuration, the flexbeam and blade properties will be different. In this thesis, 
it is assumed that the blade and the flexbeam have same uniform properties. The rea- 
sons for this assumption are bv'o fold; (i) practical blade data is not available, (ii) by 
considering uniform and same properties for both the hingeless and bearingless rotor 
blades, one can obtain a good physical insight on the effects of non-linear elastomer 
and kinematic constraint on the aeroelastic behaviour of the rotor blade. 

Even though uniform properties have been assumed for the flexbeam and blade, 
the natural frequencies of the rotating blade in flap, lag auid torsional modes are repn 
resentative of a realistic rotor blade configuration. Therefore, it is expected that this 
study will provide useful and meaningful conclusions. 

The uncoupled rotating mode shapes in flap, lag and torsion are used as the as- 
sumed modes. In determining the uncoupled rotating mode shapes, the effect of the 
elastomer is included as a linear spring element between the blade and torque tube; and 
the kinematic constraint introducing pitch-flap coupling is neglected. In the transient 
response analysis, six modes have been considered which correspond to three modes 
in flap, two modes in lag and one mode in torsion. Following the solution procedure 
outlined in Sec. 5.3, the transient response of a soft-in-plane rotor blade is obtained 
numerically incorporating a time-step of 0.003 s. These numerical results are presented 
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below in three sections. 

(i) In the first section, the blade response and the root loads are discussed for a par- 
ticular configuration treated as the baseline configuration of the rotor blade. 

(ii) The second section focusses on the influence of the locations of the elastomer and 
torque tube on the transient response characteristics of the blade. 

(iii) The third section is devoted to the effect of the location of pitch link point P 
(defined by parameters a and 6) on the transient response of the blade. When the 
parameter a is positive, the pitch link is in the leading edge side of the blade and when 
a is negative, it is in the trailing edge side of the blade. The parameter h represents the 
distance of the pitch link from the junction of blade-flexbeam-torque tube as shown in 
Fig. 5.6. 

Since the blade pitch is provided by moving the pitch control point P up or down, 
the control input chosen for the transient response is the step displacement Wp at the 
point P in the direction. 

5.4.2. 1 Baseline Configuration 

For the baseline configuration, the location of the elastomer and torque tube are set 
at = 0.10 and ^2 = 0.25 where and ^2 are shown in Fig. 5.6. The pitch link 
is considered to be located at the leading edge side of the blade. The location of the 
pitch control point P, defined by the parameters o and b are set equal to a = 0.05 and 
b = 0.15. The pitch link is given a step displacement Wp = 0.00624 in direction. 
In Figs. 5.12a-c, the tip responses of the blades in flap (wup), lag {vap) and torsional 
modes, respectively, are plotted as functions of the non-dimensional time r . The 
flap response of the baseline rotor blade reaches a steady state value 0.025 as shown 
in Fig. 5.12a. It is intersting to note that unlike the flap response of a hingeless blade 
(Fig. 5.7a), the flap response of a bearingless blade does not exhibit any.dynamic over- 
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shoot even though there exists initial oscillatory behaviour. The qualitative difference 
between the flap response of these two types of blade configurations may be attributed 
to the presence of multiple load path and kinematic constraint in the bearingless blade. 

As indicated in Fig. 5.12b, the tip response in the lag mode (vap) shows an oscilla- 
tory behaviour with a mean value of (vapj^^^n = -0.007. The amplitude of oscillation 
in the lag mode is small in comparison to that of a hingeless blade (Fig. 5.7b). The 
reduction in the response amplitude can be attributed to the presence of the elastomer 
damping in addition to the aerodj-namic damping. 

The torsional response deformation of the bearingless rotor blade at any cross- 
section is defined by the total pitch angle which is the sum of the elastic twist of the 
blade and the elastic twist at the end of the flexbeam produced by the control input 
Wp. Figure 5.12c shows the variation of the tip response in the torsional mode as 
a function of r. It is observed that the torsional response reaches a steady value 0.13 
beyond time r = 25. The presence of higher firequencies is observed in the torsional 
response. The variation of the inflow velocity (A) with r is shown in Fig. 5.13. The 
inflow velocity for bearingless blade cofiguration, as compared to that of the hingeless 
blade configuration (Fig. 5.8), takes longer time to reach the steady state value and 
no overshoot is observed. 

The variation of the constraint force Fp (acting at the point P of the pitch link) 
as a function of r is shown in Fig. 5.14. The force Fp attains the steady state value 
of 0.037 without any overshoot. The qualitative nature of the curve is same as that of 
the torsional response, shown in Fig. 5.12c. 

The loads at the flexbeam root are obtained by integrating the distributed forces 
and moments over the span of the blade. The variation of different root loads are 
shown in Figs. 5.15-5.19. It is seen from Fig. 5.15, that the out of plane root shear 
force exhibits an overshoot of 1.3 times the steady state value 0.0087. In addition. 
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it undergoes oscillation of moderate magnitude before reaching the steady state. The 
same figure also represents the response of rotor thrust coefficient (Ct) in a different 
scale. It is interesting to observe that unlike hingeless rotor, in this case, even though 
the flap response does not show any overshoot, rotor thrust exhibits overshoot (Fig. 
5.15). Figure 5.16 shows the out of plane bending moment (flap bending) of the blade. 
No appreciable overshoot is observed in this case. For the sake of completeness, the 
variation of other root loads are shown in Figs. 5.17-5.19. The in-plane root shear and 
root moment exhibit high frequency content in their response, as shown in Figs. 5.17 
and 5.18, r.espectively. The torsional moment shows a very mild oscillatory' behaviour 
as shown in Fig. 5.19. 

5A.2.2 Influence of Varying Control Input Wp 

The transient response of the blade for the baseline configuration is obtained with 
different values for the step control inputs. These values are Wp = 0.00624, 0.00936 
(1.5 X 0.00624) and 0.01248 (2.0 x 0.00624). The tip responses of the blade in flap, 
lag and torsion, corresponding to these inputs are shown in Figs. 5.20a-c. It is evident 
from the figures that irrespective of the magnitude of the input, the qualitative nature 
of the response remains same for all the modes. The steady state responses in the flap 
mode for the three input values are 0.027, 0.042 and 0.058, respectively (Fig. 5.20a). 
It is evident from the values that the steady state response in the flap exhibits a mild 
non-linear characteristics with respect to change in input values. Similar mild non- 
linear behaviour of the steady state responses in lag and torsional modes with respect 
to control pitch input Wp are observed in Figs. 5.20b and 5.20c. 
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5.4.2.3 Influence of Elastomer Location (^j, ^ 2 ) 

In this section, the influence of the elastomer and torque tube locations (^1 and ^ 2 ) on 
the transient response characteristics of the rotor blade is examined. The location of 
the control point P of the pitch link is set at a = 0.05 and 6 = 0.15. Three different 
torque tube locations (^2 = 0.20, 0.25 and 0.30) are considered. For each value of ^ 2 , 
three different elastomer locations (0.02, 0.04 and 0.10) are considered. The control 
step input w% =0.00624 is same for all the cases. 

In Fig. 5.21, the transient responses in flap, lag and torsional modes are shown 
for ^2 = 0.20 with different locations of the elastomer (^1 = 0.02, 0.04 and 0.10). The 
qualitative nature of the response in each case is found to be the same. Out of three 
configurations considered, the steady state flap and torsional responses are maximum 
for = 0.10. The amplitude of oscillation in the lag mode is also maximum for 

= 0.10. Since the distance between the elastomer and torque tube locations is least 
for this configuration, the elastomer undergoes least deformation thus providing least 
damping and stiffness to the system. Therefore, the responses exhibit higher values. 

Figures 5.22 and 5.23 present the response results corresponding to torque tube 
locations ^2 = 0.25 and ^2 = 0.30, respectively. No appreciable change is observed in the 
qualitative nature of the blade responses for these configurations; but the magnitude 
of the response is reduced with the increase in ^2 as can be seen by comparing the 
steady state values in Figs. 5.21-5.23. Comparing the lag responses for the blades with 
three different torque tube locations ^2 = 0.20, 0.25 and 0.30 (Figs. 5.21b, 5.22b and 
5.23b), it is observed that the amplitude of oscillation is less for the combinations = 
0.04, 1^2 = 0.25, and = 0.10, ^2 = 0.25 as compared to those for other configurations. 
Based on the numerical values of these results, it can be stated that the influence of 
the elastomer and torque tube locations on the transient response in the lag mode of 
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the blade is more significant compared to flap and torsional response.. 

5.4.2.4 Influence of pitch link location (a, b) 

A parametric study is conducted to examine the effect of geometric parameters a and 
b on the transient response characteristics of the blade under different pitch control 
inputs Wp. Towards this end, four diflferent v’alues of the parameter a are chosen. 
These are: a = — 0.02, 0.02, 0.03 and 0.05. The parameter b is also varied with b = 
0.10, 0.15 and 0.23. The elastomer and torque tube locations are set at = 0.10 and 
6 = 0.25. 

The transient response characteristics for the configuration a = 0.03 and b = 0.25 
are shown in Figs. 5.24a-c for three different \'alues of control input - 0.00507, 
0.0076 and 0.01014. Comparing these results \^dth those shown in Figs. 5.20a-c, it is 
observed that the transient response characteristic of the blade in all the three modes 
is similar to that for the baseline configuration (a = 0.05, b = 0.15). The results 
obtained for other combinations of a and b (not presented for conciseness) reveal that 
the qualitative nature of the blade response remains the same so long as the pitch link 
is located at the leading edge side. 

A comparison of the transient response characteristics with the pitch link at the 
leading (a = 0.02) and the trailing (a = — 0.02) edge sides is shown in Figs. 5.25-5.27. 
The parameter b is set equal to 0.15. The magnitude of the control input wf, is identical 
for the two cases. For the leading edge configurtion, the control input is a positive step 
displacement wf, = 0.00383 (e, direction) and for the trailing edge side the control 
input is set wf, = -0.00383 (-e, direction), resulting in a positive pitch input for the 
blade in both cases. 

The transient responses in the flap (wup) and lag (vap) modes are shown in Figs. 
5.25 and 5.26, respectively. FVom Fig. 5.25, it is interesting to note that the nature 
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of the flap response changes drastically depending on the sign of the value of a. The 
flap response of the blade with the trailing edge configuration exhibits an overshoot 
with a value of {wtip)max = 0.045 which is about three times the steady state value of 
{^tip)s 0.016. In contrast, no o\ershoot is observed for the leading edge configurarion. 
But the steady state values in both cases are about the same (around 0.017). 

The lag response (Fig. 5.26) at the tip, for the trailing edge configuration shows a 
large variation with a mean value (r;tip),„ean = -0.007. The variation for the leading 
edge configuration is considerably less with a mean value {vupjmean = -0.0035. 

The torsional response of the blade is represented by the response at the blade 
section at a distance 0.75L from the root. This particular blade section is usually 
chosen as the typical section for aerodynamic loads. It can be seen from Fig. 5.27 that 
the torsional response (<^ 0 . 75 ) for the trailing edge configuration exhibits an overshoot 
of about one and half times the steady state value 0.10. The torsional response for the 
leading edge configuration reaches a steady state value 0.075 with no overshoot. It is 
important to note, that for the same control input, the blade experiences a larger pitch 
deformation (<^ 0 . 75 ) for the trailing edge than that for the leading edge configuration. 

Figure 5.28 shows a substantial difference in the variation in inflow for these two 
configurations. For the trailing edge configuration, the inflow reaches a steady state 
value A = 0.05 which is higher than that of the leading edge configuration which is 
A = 0.038. 

The magnitudes of the pitch link load (| Fjp |) are plotted as a function of time (r) 
for both the configurations in Fig. 5.29. The response for the trailing edge case exhibits 
a rapid rise before reaching a steady state value (Fp), = 0.078. On the other hand, 
the pitch link load for the leading edge configuration gradually attains the steady state 
value (F», = 0.06. It may be noted that the sign of the pitch link load is negative in 
case of trailing edge pitch link configuration of the blade. 
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In order to highlight the difference between the two configurations, the root loads 
have also been evaluated. These are shown in Figs. 5.30-5.34. The variation of the 
out of plane root shear as a function of time is shown in Fig. 5.30. Both the blade 
configurations exhibit overshoot. But for the trailing edge configuration, the overshoot 
is very large, i.e., about five times the steady state value 0.006. For the leading edge 
configuration, the overshoot is about one and half times the steady state value 0.004. 
This steady state value of the out of plane root shear for the trailing edge coonfiguration 
is 1.5 times that for the leading edge configuration. This increase is due to the effect 
of increased torsional response and inflow velocity as shown in Figs. 5.27 and 5.28, 
respectively. The out of plane root moment shown in Fig. 5.31, indicates that the 
overshoot for the trailing edge configuration is about four times the steady state value, 
while the overshoot for the leading edge is very small (about 1.1 times). 

The variation of the in-plane root shear and moment are shown in Figs. 5.32 and 
5.33, respectively. For the trailing edge configuration, the response seems to contain 
high frequency components. The amplitudes of the loads is much larger for the trailing 
edge configuration as compared to that for the leading edge configuration. 

The torsional root moments for both configurations are shown in Fig. 5.34. There 
is a rapid rise in the torsional moment for the trailing edge at the initial stage whereas 
the rise is gradual for the leading edge configuration. The nature of the response is 
very similar to the torsional response of the blade as shown in Fig. 5.27. Based on 
these results it can be stated that there is a qualitative change in the nature of the 
response depending on whether the pitch link is at the leading or trailing edge. 

It had been observed that for certain trailing edge configurations, the blade response 
did not converge. To highlight this behaviour, the blade response for the two cases 
corresponding to (i) leading edge configuration with a = 0.03 and b = 0.23, and (ii) 
trailing edge configuration with a = - 0.03 and b = 0.23, have ben evalpted. Control 
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input Wp = 0.00525 and Wp = —0.00525 are applied for the leading and trailing edge 
configurations, respectively. The transient response in the flap, lag and torsional modes 
are shown in Fig. 5.35 for both these configurations. The blade response reaches 
a steady state value for the leading edge pitch link configuration, but the response 
diverges for the trailing edge configuration indicating instability. 


5.4.2.5 Aeroelastic Coupling 


Aeroelastic couplings have a strong influence on the response and stability character- 
istics of a rotor blade. In the case of kinematic pitch-flap and pitch-lag couplings, 
quantitative measure of the couplings are defined based on the orientation of the lag 
or flap hinges [113]. In hingeless and bearingless rotors, due to the absence of hinges. 
It is difficult to quantify how much of blade torsional motion is due to elastic twist 
and how much is due to aeroelastic pitch-flap and pitch-lag couplings. Therefore the 
quantitative measure of the couplings have to be defined based on the blade response. 
Since the response depends on the blade loading, the couplings will also be dependent 
on the blade loading. In this section, using the steady state tip response of the blade, 
an attempt is made to identify a measure of the aeroelastic coupling present in the 
bearingless rotor blade system. 

In Ref. [48] Hodges and Ormiston proposed a mathematical expression to measure 
the equivalent kinematic pitch-flap and pitch-lag couplings generated in elastic blades. 


Following this idea, the pitch-flap (kf}) and pitch-lag coupling (A:^) are defined as 


k^ = 


and k( = 



(5.70) 

(5.71) 


where subscript s denotes the steady state response and 


00.75 = the steady state value of the torsional response at 0.75L of the blade. 
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Wttp — the steady state tip response in the flap mode, 


and vtip = the steady state response in the lag mode. 

Since the convergence in lag mode is very slow, the steady state value {vtip)s is taken 
as the mean value of the oscillatory lag response. This mean value is computed from 
the oscillatory lag response of the blade. The influence of the geometric parameters a, 
h and the magnitude of the control input on the aeroelastic coupling measures, defined 
by Eqs. 5.70 and 5.71, are evaluated. The step input of the control link displacement 
Wp is varied within a range of w], = -0.00255 to 0.01248. The elastomer location is 
set at = 0.10 and ^2 = 0.25. The values of the parameter a are chosen as a = - 
0.02, 0.02, 0.03 and 0.05. In addition, three different values of the parameter b are also 
considered: b = 0.10, 0.15 and 0.23. 

The steady state response of the blade are obtained for several combinations of 
the parameters. In each case, the aeroelastic couplings are calculated using the steady 
state response of the blade. The numerical values of the couplings as a function of 
control input and parameter b are presented in Tables 5.2 and 5.3. Table 5.2 presents 
the results corresponding to the case a = 0.03 and Table 5.3 provides the results for 
the case a = 0.05. It is seen from Table 5.2 that with increase in the control input 
w’p, the magnitudes of the pitch-flap (kp) and pitch-lag {k^) couplings decrease. The 
reduction in the magnitude is more prominent in the case of pitch-lag coupling than for 
pitch-flap coupling. A similar behaviour is also observed for the case a = 0.05 as shown 
in Table 5.3. A pictorial representation of the variation of the aeroelastic couplings are 
also presented in 3-D plots. (Figures 5.36 and 5.37 correspond to the case a = 0.3 and 

Figs. 5.38 and 5.39 pertaining to the case a = 0.05.) 

A comparison of the aeroelastic coupling measures between the leading and the 

trailing edge configurations is presented in Table 5.4. The two configurations (i) 
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a = - 0.02, b = 0.15, and (ii) a = +0.02, 6 = 0.15. The aeroelastic couplings are 
evaluated for three values of the control inputs w^p. It is observed from Table 5.4 that 
for both configurations the magnitudes of pitch-flap and pitch-lag couplings decrease 
with increasing control input. For both configurations, the reduction in pitch-flap 
coupling is found to be more pronounced than the pitch-lag coupling. In addition, 
the reduction in coupling values seems to be more for the trailing edge configuration 
than the leading edge configuration. While the pitch-flap coupling for trailing edge 
configuration is higher than the leading edge configuration, the magnitude of the pitch- 
lag coupling is lower for trailing edge configuration than the leading edge configuration. 
The reason for the difference in aeroelastic couplings for these two configurations may 
be attributed to the change in the direction of the constraint force and moments with 
a change in the sign of parameter a given in Eq. 5.54 and 5.58. 

The dependence of pitch-flap and pitch-lag coupling with control input can be at- 
tributed to the non-linear elastomer and other non-linear eflfects in the blade dynamics. 
The results clearly indicates that the aeroelastic coupling in bearingless main rotor is 
dependent on the operating condition of the blade and also on the other geometric 
configuration. 

5.5 Summary 

A suitable numerical technique has been formulated to study the transient response of 
an isolated bearingless rotor blade to a step control pitch input. The blade is assumed to 
undergo coupled flap-lag-torsional deformation in hovering condition. A time varying 
inflow model (based on the dynamic inflow model) has been used for the analysis. 
The expressions for aerodynamic loads are not obtained explicitly, rather they are 
derived an implicit formulation. The study is carried out for different locations of the 
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elastomer and pitch link. The transient response results of the bearingless rotor blade 
are compared with those of the hingeless blade. The comparison is made to highlight 
the difference in their response behaviour. It is noticed that the transient response 
characteristics of a bearingless blade depends on the pitch link locations. In particular, 
it is observed that the blade response with trailing edge pitch link configuration is 
significantly different from that of a leading edge configuration. For certain trailing edge 
configurations, the blade exhibits a divergence transient response indicating instability. 
A quantitative measure of the aeroelastic couplings (pitch-flap, pitch-lag), based on 
the steady state response has been evaluated and their dependence on the geometric 
configurations and control pitch input has been identified. 
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Fig. 5.2 Undeformed cross-section fixed rotating coordinate system. 
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Fig. 5.3 Undeformed curvilinear cross-sectional coordinate system. 



Fig. 5.4 Deformed cross-sectional coordinate system and components of blade 
velocity. 
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Fig. 5.5 Components of sectional aerodynamic forces. 
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Fig. 5.6 Bearingless rotor blade with pitch link at the leading edge side. 
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Fig. 5.10 Variation of CTmaxlCT$ for hingeless blade for various rates of pitch input, 
o different pitch rates. 
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Fig. 5.11 Comparison of experimental and calculated ratios of CTmaxICrs 
hingeless blade for various rates of pitch input. 

calculated curve; □, o A, Experimental data. (Taken from Ref. [106]). 
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Fig. 5.13 Variation of inflow for bearingless blade to a step control pitch input 
(wj, = 0.00624). 

. Configuration: = 0.10, ^2 = 0.25 and a = 0.05, b = 0.15. 


Fig. 5.14 Variation of constraint force Fp for bearingless blade to a step control pitch 
input {wp = 0.00624). 

Configuration: = 0.10, ^2 = 0.25 and a = 0.05, b = 0.15. 
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Fig. 5.15 Variation of out of plane root shear and thrust for bearingless blade to a 
step control pitch input (wf, = 0.00624). 

Configuration: = 0.10, ^2 = 0.25 and a = 0.05, b = 0.15. 




Fig. 5.16 Variation of out of plane root moment for bearingless blade to a step control 
pitch input {wp = 0.00624). 

Configuration: = 0.10, ^2 = 0.25 and a = 0.05, b = 0.15. 
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Fig. 5.18 Variation of in-plane root moment for bearingless blade to a step control 


pitch input {wp = 0.00624). 

Configuration: = 0.10, ^2 = 0-25 and a = 0.05, b = 0.15. 


149 





r 


Fig. 5.19 Variation of torsional root moment for bearingless blade to a step control 
pitch input {wp = 0.00624). 

Configuration: = 0.10, ^2 = 0.25 and a = 0.05, b = 0.15. 



150 




0.07 

0.06 

0.05 


25 30 35 40 


Umi 


0 5 10 15 20 25 30 35 40 


0.05 h 



(c) Torsion 
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Fig. 5.20 Bearingless blade response with different values for step control pitch 
inputs. w‘p = 0.00624; = 0.00936; w^p = 0.01248. 

Configuration; = 0.10, ^2 = 0-25 and a = 0.05, b = 0.15. 









Fig. 5.21 Bearingless blade response with different locations of the elastomer (^i) 
6=0.04; 6=0.02; — 6 = 0.1. 


Configuration: 6 = 0.2, Input: Wp = 0.00624; a = 0.05, b = 0.15. 






-e- 



(c) Torsion 


30 35 40 


Fig. 5.22 Bearingless blade response with different locations of the elastomer (^i). 
ei=0.04; 6=0.02; 6=0-1. 


Configuration: 6 = 0.25, Input; Wp = 0.00624; a = 0.05, b = 0.15. 
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(c) Torsion 



Fig. 5.23 Bearingless blade response with different locations of the elastomer (^i). 

6 = 0.1; 6 = 0.02; — 6 = 0.4. 

Configuration: ^2 = 0.3, Input: Wp = 0.00624; a = 0.05, b — 0.15. 






(c) Torsion 
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Fig. 5.24 Bearingless blade response with different values for step control pitch 

inputs. wj> = 0.00507; w‘p = 0.00761; w^p = 0.01014. 

Configuration: = 0.10, ^2 = 0.25 and a = 0.03, b = 0.15. 






Fig. 5.25 Bearingless blade response in flap mode with different pitch link locations. 


a = -0.02 and = -0.00383; a = 0.02 and w^p = 0.00383. 

Configuration: = 0.10, ^2 = 0.25 and b = 0.15. 



Fig, 5.26 Bearingless blade response in lag mode with different pitch link locations. 

a = —0.02 and Wp = -0.00383; — a = 0.02 and Wp = 0.00383. 

Configuration: ^1 = 0.10, ^2 = 0.25 and h = 0.15. 
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Fig. 5.27 Bearingless blade response in torsional mode with different pitch link 

locations. a = -0.02 and w^p = -0.00383; a = 0.02 and w‘p = 0.00383 

Configuration: = 0.10, ^2 = 0.25 and b = 0.15. 



r 

Fig. 5.28 Variation of inflow for bearingless blade with different pitch link locations. 

0 = -0.02 and w^p = -0.00383; a = 0.02 and wj, = 0.00383. 

Configuration: = 0.10, ^2 = 0.25 and b = 0.15. 






T 


Fig. 5.29 Variation of constraint force Fp for bearingless blade with dilferent pitch 
link locations. 

a = -0.02 and Wp = —0.00383; a = 0.02 and Wp = 0.00383. 

Configuration; = 0.10, ^2 = 0-25 and h = 0.15. 
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Fig. 5.30 Variation of out of plane root shear for bearingless blade with different pitch 
link locations. 

a = -0.02 and wj, = -0.00383; a = 0.02 and w}, = 0.00383. 

Configuration: = 0.10, ^2 = 0.25 and b = 0.15. 
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Fig. 5.31 Variation of out of plane root moment for bearingless blade with different 
pitch link locations. 

a= -0.02 and wj> = -0.00383; a = 0.02 and w‘p = 0.00383. 

Configuration: = 0.10, ^2 = 0.25 and b = 0.15. 
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Fig. 5.33 Variation of in-plane root moment for bearingless blade with different pitch 
link locations. 

a = —0.02 and Wp = —0.00383; a = 0.02 and Wp = 0.00383. 

Configuration: = 0.10, ^2 = 0.25 and b = 0.15. 
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Fig. 5.34 \ ariation of torsional root moment for bearingless blade with different pitch 
link locations. 

a = -0.02 and = -0.00383; a = 0.02 and wf> = 0.00383. 

Configuration: = 0.10, ^2 = 0.25 and b = 0.15. 
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Fig. 5.36 Influence of pitch link location (6) and step control pitch input {wp) on 
pitch-flap coupling of a bearingless blade. 

Configuration: = 0.10, ^2 = 0.25 and a — 0.03. 



Fig. 5.37 Influence of pitch link location (b) and step control pitch input (wp) on 
pitch-lag coupling of a bearingless blade. 

Configuration: = 0.10, ^2 = 0.25 and a = 0.03. 



Fig. 5.38 Influence of pitch link location (b) and step control pitch input {wp) on 
pitch-flap coupling of a bearingless blade. 

Configuration: = 0.10, ^2 = 0-25 and a = 0.05. 



Fig. 5.39 Influence of pitch link location (6) and step control pitch input (w’p) on 
pitch-lag coupling of a bearingless blade. 

Configuration: = 0.10, ^2 = 0.25 and a = 0.05. 
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Table 5.1 Input data for transient response analysis. 


7 


5 

k 


TT 

L 


20 

Cl 


2it 

Cdo 


0.01 

N 


4 

GJ 

mCl^L* 


0.001473 



0.0301 

SIfjrj 

mQ^L^ 


0.0106 

Inirj^ 

UnJJ 


0.0004 



0.0 



0.0 

(tr 


1.5 

k\ = ^ 

EA ’ - mL^ 



Offsets of mass centre, tension centre and aerodynamic centre from the elastic axis are 
zero. 
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Table 5.2 Variation of aeroelastic coupling with b and step control pitch input Wp. 
= 0.1, ^2 = 0-25 and a = 0.03. 


b 

Wp 


k 

0.10 

0.00361 

4.39 

22.17 

0.10 

0.00542 

4.02 

20.67 

0.10 

0.00722 

3.84 

19.86 

0.15 

0.00507 

4.50 

18.98 

0.15 

0.00761 

4.08 

15.59 

0.15 

0.01014 

3.80 

13.03 

0.23 

0.00525 

5.04 

14.24 

0.23 

0.00787 

4.45 

11.31 

0.23 

0.01050 

4.05 

9.28 




Table 5.3 Variation of aeroelastic coupling with b and step control pitch input Wp. 
= 0.1, ^2 = 0-25 and a = 0.05. 


b 

Wp 

kp 


0.10 

0.00537 

4.40 

20.41 

0.10 

0.00805 

4.01 

17.53 

0.10 

0.01070 

3.74 

14.71 

0.15 

0.00624 

4.48 

16.58 

0.15 

0.00936 

4.03 

13.73 

0.15 

0.01248 

3.74 

11.01 

0.23 

0.00475 

5.15 

15.40 

0.23 

0.00712 

4.57 

11.88 

0.23 

0.00950 

4.22 

10.04 




Table 5.4 Variation of aeroelastic coupling with a and step control pitch input Wp. 
= 0.1, ^2 = 0-25 and b = 0.15. 


a 

Wp 


k 

-0.02 

-0.00255 

7.79 

-16.59 

-0.02 

-0.00383 

6.12 

-15.98 

-0.02 

-0.00511 

5.88 

-15.07 

0.02 

0.00255 

5.01 

-20.28 

0.02 

0.00383 

4.63 

-20.21 

0.02 

0.00511 

4.35 

-20.01 
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Chapter 6 


ELASTOMERIC DAMPER 
MODEL AND LIMIT CYCLE 
OSCILLATION 

6.1 Introduction 

The review of the literature given in Sec. 1.2 indicates that some of the recent studies 
have focussed on the phenomenon of limit cycle behaviour of the rotor blades due to the 
elastomer nonlinearity. Ormiston et al. [31] considered an elastomer model in which 
the damping force was assumed to be proportional to a linear combination of different 
powers of velocity (powers of 1/2, 1, 2 and 3). Using this non-linear elastomer model, 
they showed the presence of limit cycle oscillation for both hingeless and bearingless 
rotor blade configurations. It was also pointed out that there was no limit cycle when 
the elastomer damping is taken as linear. 

The hover air-resonance test conducted by Panda et al. [4] also indicated that 
non-linearities in stiffness and damping characterisrics of the elastomer seem to have 
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significant influence on the limit cycle oscillations. When they used a fluidlastic damper 
having linear characteristics, no limit cycle oscillation was observed. In Ref. [30], 
Gandhi and Chopra proposed a non-linear visco-elastic model for the elastomer in 
which the damping element is linear. The model shows a rapid decrease in the damping 
coefficient at very low amplitudes of oscillation; but its value is always positive. With 
this elastomeric damper, they showed a limit cycle oscillation for an autonomous system 
representing the isolated lag dynamics of a blade. 

The theorem on the existence of a limit cycle [121] says that for the equation 

x + f {x)x + g{x) == 0 (6.1) 

to exhibit a limit cycle, f{x) must be positive when | x | is large, and negative when 
I 2 : 1 is small, and where g is such that, in the absence of the damping term f{x)x, a 
periodic solution is expected for small x. In Van der Pol’s equation 

X -I- e - l) i -1- a: = 0, e > 0 (6.2) 

j.e., the damping coefficient becomes negative in the strip | x |< 1 and positve for 
I X |> 1. Hence the system can show a stable limit cycle for any given initial distur- 
bance. Therefore, it is difficult to physically understand as to why a damper whose 
coefficient is always positive would lead to a limit cycle oscillation under free vibration 
conditions as reported in Ref. [30]. For such a strictly dissipative system, the energy 
source for sustained oscillation is absent. To explore the phenomenon of limit cycle, 
in this chapter a study is carried out on the response of an isolated bearingless rotor 
blade using different elastomer models. 
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6.2 Formulation of Elastomer Model 


Since experimental data are used in developing the analytical spring-damper model for 
the elastomer, there is a scope for formulating many different types of models. In Chap- 
ter 2, the elastomer is modelled (herein denoted as Model 1) as a parallel combination 
of a non-linear spring, a Coulomb damper and a hysteretic damping. A modified model 
(denoted as Model 2) consisting of a non-linear spring, a Coulomb damping element 
and a Rayleigh type hysteretic damping element can also be considered (Fig. 6.1). The 
constitutive differential equation of the elastomer model under harmonic loading for 
Model 1 and Model 2 can then be written, respectively, as 

Kix — Kzx^ -f K^x^ — Kjx^ + Fsgn | x | -l — -x = DoSincji Model 1 (6.3) 

u 

KiX — K^x^ -I- K^x^ — Kjx^ + Fsgn | x | -I — ^x^ — —x = Do sin wt Model 2 (6.4) 

Assuming a harmonic motion and following the procedure mentioned in Sec. 2.2, 
the in-phase stiffness {G') and the quadrature stiffness {G") for Model 2 can be obtained 
as 


G' = 

ATi - - IkjX^ 

4 8 2 

(6.5) 

and G" = 

+ 7^3-^ ~ • 

TTA 4 

(6.6) 


6.3 Results and Discussion 

Solving the minimisation problem described in Chapter 2, the parameters for the elas- 
tomer Model 2 are obtained. These are given in Table 6.1 along with the parameters 
of Model 1. Comparison of the experimental data with Model 1 and Model 2 is shown 
in Fig. 6.2. It can be seen that both the models fit the experimental values fairly well 
in the available range of data. 
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Using the Model 1 and Model 2, the transient response characteristics of the rotor 
blade (coupled flap, lag and torsion) are obtained following the solution procedure 
outlined in Sec. 5.3. The blade data is given in Table 5.1. The elastomer and torque 
tube locations are set at = 0.10 and ^2 = 0.25. The pitch link is located at the 
leading edge side of the blade with a = 0.03 and b = 0.10. A step control input of Wp 
= 0.00361 is applied. Figures 6.3 and 6.4 show the transient response in lag mode for 
the two cases. In the case of Model 2 (Fig. 6.3), it is observed that within a short 
time, the lag response settles down to a stable limit cycle oscillation with a constant 
amplitude of 0.0145. In the case of Model 1 (Fig. 6.4), the transient response in lag 
mode settles to steady state value of - 0.0027. The blade response shovm in Figs. 6.3 
and 6.4 are obtained with the inclusion of aerodynamic loading. 

In order to ensure that the cause of limit cycle is due to the elastomer model and 
not due to aerodynamic loading, the transient analysis is also carried out by excluding 
the aerodynamic effects. For the same set of blade parameters and input, the response 
of the blade is computed. Figures 6.5 and 6.6 show the responses in the lag mode of 
the blade with elastomer Model 2 and Model 1, respectively. While the response of 
Model 2 exhibits a limit cycle oscillation (Fig. 6.5), no such oscillation is observed for 
Model 1 (Fig. 6.6). 

6.4 Summary 

Even though both elastomer models fit the experimental data fairly well, they give 
rise to totally different response characteristics. These results clearly indicate that 
the phenomenon of limit cycle oscillation is highly dependent on the elastomer model. 
Hence, care must be excercised in modelling the elastomer. 


174 



X 10® (N/m) 


X 



Fig. 6.1 Elastomer Model 2 with Rayleigh type damping. 



X X 10-^ 


Fig. 6.2 Comparison of the experimental data with Model 1 and Model 2. 
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Table 6.1 System parameters of the elastomeric damper Model 1 and Model 2. 


System 

Parameter 

Model 1 

Model 2 

K, (N/m) 
Kz (N/m^) 
K, (N/m^) 
K, (N/m^) 
F(N) 
h (N/m) 
hi (N/m) 
hz (N/m3) 

2.673989 x 10® 
1.315287 X 10^2 
3.519586 X 10^^ 
3.176266 X 10^2 
4.797347 x 10^ 
4.569120 X 10® 

2.673989 x 10® 
1.315287 X 10^2 
3.519586 X 10^2 
3.176266 X 1022 
8.45909 X 102 

1.44619 X 10® 
7.06419 X 10*° 







Chapter 7 


CONCLUSIONS 


The effect of nonlinearity of the elastomer on the dynamics of an idealised bearingless 
rotor blade has been studied. Attempts have been made to analyse problems of in- 
creasing complexities with an aim to gain a fundamental understanding of the blade 
dynamics and aeroelastic transient response. In the absence practical blade data, uni- 
form properties reported in research publications are used for the blade. The properties 
of the elastomer are taken from the experimental results published in the open litera- 
ture. Even though the blade is assumed to be uniform, the natural frequencies of the 
blade in flap, lag and torsional modes are representative of an actual blade configura- 
tion. The major conclusions of the present study are summarised below. The results 
are valid within the framework of the assumptions on which the study is carried out. 

1. The experimental data indicate that the stiflFness and damping properties of the 
elastomer exhibit very weak dependence on frequency in the range of interest for rotor 
blade analysis. In view of this observation, a simple model for the non-linear elastomer 
has been developed which can be easily integrated with the blade model. The highly 
non-linear characteristic of the elastomer is modelled by a parallel combination of spring 
and damper elements. 
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2. Using an idealised model for a beaingless rotor, both linear and non-linear free 
vibration analyses of the rotating blade under uncoupled flap and lag bending have 
been studied. For the linear problem, two different solution techniques, one based on 
power series expansion and the other based on the Rayleigh-Ritz method are used. 
The natural frequencies obtained from these two methods are found to be in excellent 
agreement. For the flap motion, it is observed that the effect of the elastomer spring 
has a relatively high influence on the third mode; a moderate influence on the first 
mode; and less influence on the second mode. In case of lag motion, the eflfect of the 
spring is more pronounced on the first mode. 

3. In the non-linear analysis, a numerical-perturbation technique based on multiple- 
time-scale is formulated to determine a closed form frequency-amplitude relationship 
for uncoupled flap and lag motions of the rotating blade. It is concluded that upto a 
fairly high value of amplitude, not normally to be exceeded in practice, a seventh order 
expansion of the non-linear spring force is sufficient to correctly predict the frequency- 
amplitude relationship. If the amplitude level is lower, one can use an expansion upto 
fifth order and still get accurate results while saving on the computational effort. The 
alternating signs in the polynomial, expressing the restoring force prevents the non- 
linear characteristic to be monotonically softening or hardening. 

4. An amplitude dependent stability anlysis has been j^rformed for a coupled 
rotor /fuselage system under ground resonance condition. It is observed that the elas- 
tomer amplitude does not significantly alter the stability of the system at the region of 
maximum instability. But the damping in the progressive lag mode decreases appre- 
ciably with the increase in amplitude. 

5. The locations of the elastomer and torque tube seem to have significant influence 
on the stability of the system. The results indicate that for a given location of the 
elastomer, there is an optimum location for the attachment of the torque tube providing 
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least instability in the system. The physical reason for this observ^ation is that for a 
particular combination of elastomer and torque tube locations, the deformation of 
the elastomer is maximum, resulting in a maximum value of damping. Therefore, a 
judicial selection of the locations of the elastomer and torque tube is very important 
for enhancing the stability of the system under ground resonance condition. 

6. A suitable numerical technique has been formulated to study the transient 
response, to a step input, of an isolated bearingless rotor blade undergoing coupled flap- 
lag-torsional deformation in hovering condition. A similar analysis is also performed 
for a hingeless blade to validate the equations of motion and solution technique. It 
is interesting to note that there is a qualitative difference in the nature of the blade 
response and rotor thrust for these two types of blades. The reason for the difference 
may be attributed to the presence of multiple load path and kinematic constraint in 
the bearingless blade. 

7. It is also observed that the transient response characteristics change drastically 
depending on whether the pitch link is located at the leading edge or the trailing edge 
side of the blade. The root loads for these configurations also exhibit similar differences. 
These differences may be attributed to the change in the directions of the constraint 
force and moments when the pitch link is shifted from leading to the trailing edge side. 

8. An attempt is made to define a measure of the aeroelastic coupling present 
in the bearingless rotor blade system. The influence of the geometric parameters a, 
b and the magnitudes of the control input on the aeroelastic coupling measures are 
evaluated. The magnitude of the pitch-lag coupling is found to be higher for the 
leading edge configuration than that for the trailing edge configuration. The trend is 
found to be reversed for the pitch-flap coupling. The aeroelastic coupling measures are 
observed to be dependent on the magnitude of the control input, exhibiting a non-linear 
relationship. 
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9. A study is carried out on the response of an isolated bearingless rotor blade using 
two different elastomer models representing the same experimental data. The results 
indicate that the phenomenon of limit cjmle oscillation is highly dependent on the 
elastomer model. Hence, sufficient care must be exercised in modelling the elastomer. 
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Appendix A 


Derivation of Strain Expressions 

The procedure for deriving the strain expressions mentioned in Sec. 5. 2.4.1 follows the 
method outlined in Ref. [46]. 

A.l Coordinate Transformation 

For the derivation of strain expressions, it is essential to establish transformation rela- 
tions between various coordinate systems used in this formulation. 

The transformation relation between the undeformed curvilinear cross-sectional 
coordinate system (ex,e,„e^) and the cross-section fixed rotating coordinate system 
{exi,eyi,ezi) can be defined as (Fig. 5.3) 

(A.l) 




6x 



. = [Tx]< 


6x1 

Cyl 

ezi 
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where 


1 0 0 

= I 0 cos 6 sin ^ I (A. 2) 

0 — sin 6 cos 6 

Since the system (ixi,^yi,^zi) is parallel to (ci, e^, iz) system, the transformation given 
by Eq. A. 2 is also applicable between (4,e,„e^) and (4,ey,6^) 

Differentiation of Eq. A.l with respect to x yields 


< > 


0 


► = 

Toe^ ‘ 





0 

0 

0 


Cx 

0 

0 

To 

i 

ir, 

0 

-To 

0 




where tq = 6^x 


(A.3) 

(A.4) 


represents the initial twist of the blade. 

The transformation between the deformed curvilinear coordinate system and the 
undeformed curvilinear coordinate system is due to rigid body rotation. It can be 
defined as 


f '' 

€'x^ 


ex 

erf 

■ = m . 






where the transformation matrix [T 2 ] is given by 


(A.5) 



1 

0 

0 


cos 6,, 

0 

sin Oj, 


COS Oq 

sin Or, 

0 

\T2] = 

0 

cos 6x 

sin Ox 


0 

1 

0 


— sin Or, 

cos 0(^ 

0 


0 

— sin dx 

cos Ox 


- sin Or, 

0 

cos Or, 


0 

0 

1 


(A.6) 
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The above transformation is obtained by following lag-flap-torsion sequence of motions. 
First the blade is rotated by about . Then it is rotated by about the rotated 
axis. Finally, it is rotated by 6x about the rotated Cx axis. The deformation sequence 
is shown in Fig. A.l. 

Using Eqs. A.l and A. 6, the relation between the deformed curvilinear coordinate 
system and cross-sectional fixed rotating coordinate system can be written as 


r > 

Cx* 


irf 

• = m| 

k' 

V, ^ 

[ e.i _ 


where 


(A.7) 


[T] = [r.][rj 


(A.8) 


The relationships between the Euler angles {Oj^, 6^, 4>) and the displacement variables 
(u, V, w) can be obtained from Fig. A.l as 


<i> 


sin Or, 


cos Or, = 


sin 6(^ = 


cos 6„ 


v^x sin 6 — w^x cos 0 


V^(l + u + (», X cos 0 -f- w^x sin Oy + {wx cos 0 - sin Oy 

\/(l + '^,xY + (n,! cos 6 -f Wx sin Oy 

X cos 6 + Wx sin Oy (wx cos 0 — Vx sin Oy 
(vx cos 0 + Wx sin 6 
^/(l + u^) + (u I cos 0 + Wx sin Oy 

1 + Ux 


(A.9) 


^(1+uy + (ti,, COS 6 + Wx sin Oy 
Applying ordering scheme, the transformation relationship becomes 


(T) = 




w 




—Vx cos {6 + (j)) - Wx sin {6 + 4>) cos {0 -I- <l>) sin {6 -I- 

Vx sin {9 + (j)) — Wx cos {6 + <l>) — sin (0 -t- ^) + cos 6 cos {9 + 4>) + sin 9 

(A.IO) 
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where 


t'c = {v^x sin 6 — w^x cos 9){v^x cos 6 + w^x sin 9) 
With the small angle approximation, Eq. A. 10 reduces to 


(A.ll) 


[T] = 


1 W^x 


(A.12) 


-v^x - w^x{9 + (f>) 1 {9 + (l>) 

v,x{^ + ^ - '^,xW,x) 1 

The relationship between the undeformed curvilinear cross-sectional coordinate sys- 
tem and the hub fixed rotating coordinate system is also given by 


(A.13) 


A. 2 Kinematics of Deformation 

The position vector of a point P on the undeformed beam with respect to the hub 
centre is given by 


ex' 


> 

&r,' 

>=m< 


h' 

V ’ 




rix, T], C) = x^x + -I- Cec 

The undeformed base vectors at P are 


(A.14) 


9t] ~ ~ 

9q = = 

where tq = 


(A.15) 
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The position vector of the point P in the deformed configuration can be written as 

R C) = Ro {x) + -nEr, + [t], C) Ci' (A. 16) 

where Ro{x) = R {x, 0, 0) 

= {x + u) ix + vcy ■+ wcz (A. 17) 

The first three terms of the right hand side of Ekj. A. 16 represent translations and ro- 
tations of cross-section. The last term is the out of plane warping of the cross-section 
where a is the amplitude of warping. Since the coordinate system is ob- 

tained by translation and rotation of the coordinate system (ex,e^, ej, the unit vector 
e^> can be assumed to be in the direction of Ex i.e., along the tangent to the deformed 
elastic axis of the beam, while the orientations of e^' and e^' are approximately equal 
to E,j and E^, respectively [122]. Following the definition given in Ref [122], the 
deformed base vectors at a point on the elastic axis can be expressed as 

Ejf — 2 €xti^x' 'b “b C»7>j) "b (A. 18) 

E^ = 2€xq^x' + + ^cc) ®c' 

The assumption of no sheeir deformation and no deformation of the cross-section leads 
to 

(■Tfri ~ ~ ^ 

The deformed base vector at the elastic axis can be expressed as 

Ex — (Id" Cxi) ^i' 

E, = erf (A.19) 

Eq = k' 
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The deformed base vectors at the point P are: 


Gx = = Ex + rjCrf^^ + Cec'.x + a, x^e^.x 

G, = = (A.20) 

= ■R.c = h' + 


The component of strain tensor can be defined in the curvilinear cross-sectional 
coordinate system as 

fi, = l{{G,-Gj)-{ 9 .- 9 j)} , hj = x,r),(: (A.21) 

Combining Eqs. A.15, A.19 and A.20, the strain expressions can be rewritten as 

/.X = e,.-r]k,-Ch + a.'^ + l{v^ + e){r^-r!) 

Ur, = U = \W^,r,-C{r-To)] 

Uc = /cx = ^[a^,f + ^(^-^‘>)] 
fm, « 0 
fa ~ 0 

fnc. = /<»? ~ 0 


where the derivatives of triad (lx, e,- , e^') are related to the curvatures A:,, and twist 
r of the deformed elastic axis of the beam by [45] 



Strain Component in Local Cartesian System 

The stress-strain relation of a material is generally given in the local cartesian 
system (y..».V3) -‘h uni, vectors e, U = 1,2,3). Therefore, i, is essential to express 
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the strain components in the local coordinate system. The expression for the strain 
tensor in the local cartesian coordinates can be written as [45, ] 


3 3 


idy. 

where Xj (i = 1,2,3) refer to the curvilinear 
relation between curvilin 
straight blade is given by 


^ 

coordinates system. The transformation 
relation between curvilinear system (x,)and local cartesian system (i/t) 


dxi 


By, 


= [9k - 9. 


<,]'"[9k-^i\ 


f 


-7]ro 


0 0 
1 0 
0 1 


(A.25) 


•where 


[9k-9i] = 


l + + 'Cto nro 

1 


r]TO 


0 


0 

1 


(A.26) 


and 


bfc-ej] = 


1 -(To riTo 
0 1 ^ 
0 0 


1 


(A.27) 


(ci, C 2 ) ^ 3 ) parallel to (ix, e,,, e^).) 


^,’s are given in Eq. A. 15. (It may be noted that , 

The determinant of the matric tensor [p* • 9^ represented by sy 

g = dei [Ok • 9i] 


(A.28) 


Combining Eqs. A.22, A.25 and A.24, one obtains the strain components 


local cartesian coordinates as 

JL 

Cx* - Vk,, - (k + + 2 


^XX 


i (>)" + ci 


189 



= Q'if,,} -((t- To) 
'yxc = ^ ('T - To) 


(A.29) 


- ^C< = 7t?c = 0 

where 


Txfj — 26 j;^ , yx^ — 26x^ , 7i7C — 

For the aeroelastic analysis, it is essential to express the above strain components in 
terms of the displacement components (w, v, w) of the elastic axis and the elastic twist 
angle ((^). By differentiating Eq. A. 17 with respect to x, the deformed base vector at 
the elastic axis is given by 

= (1 + W,x) Cx + V^x^y + (A.30) 


Combining Eqs. A. 19 and A.30, the magnitude of Ex can be written as 


1^ = (1 + u)^ + Ux^ + W^x^ — (1 + €xi)^ 


(A.31) 


Neglecting the terms and ('u.x)^) the axial strain €xx can be expressed in terms 
of the displacement components (n, n, w) as 


^XX 


Ux + + 2^^ 


(A.32) 


Combining Eqs. A.3, A. 6, A.9 and A. 23 and applying the consistent ordering scheme, 
expression for the deformed curvatures and twist can be given by 


krj = v,xi cos (^ + ^) + w^x sin {6 4- 4) 
h — ~‘^,xx sin {6 + 4)-^ ty,ri cos {6 + 4) 
T = Tf,-^ 4, x + 4o 


(A.33) 
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where (?i>o, the second order elastic twist effect of the beam can be given by 

(t>o = sin 9 + w^x cos 9) {v^xx cos 9 + w^xx sin 9) (A.34) 

Introducing Eqs. A. 32 and A. 33 into Eq. A. 29 and assuming the warping amplitude 
to be equal to derivative of elastic twist 

a = -<f>^x (A. 35) 

and neglecting the second order twist effect, the strain components can be expressed 
as 

Cix = Ux + -f - v^x [t] cos (^ + ^) - C sin {9 + <l>)] 

i2 

- U'^T[»7Sin(tf + ^) + Ccos(^ + ^)]+ 

1„ = (A-36) 

7ic = + ‘n<l>,x 
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